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Introduction 


Munitions  use  on  military  lands  is  of  keen  environmental  interest  primarily  due  to 
the  potential  for  contamination  from  unexploded  ordnance  (UXO)  constituents.  There 
has  been  little  research  on  ecosystem  response  to  munitions  use.  Demarais,  et  al., 
(1999)  describe  a  wide  array  of  spatial  and  temporal  disturbance  regimes  and  agents  of 
disturbance  associated  with  military  training.  They  summarized  the  limited  body  of 
knowledge  related  to  the  live-fire  training  and  its  physical,  chemical  and  biological 
effects  on  the  environment.  However,  little  is  known  of  the  ability  or  capacity  of  an 
impact  area  ecosystem  to  response  to  munitions  disturbance.  Munitions  use  in  the 
Washington  Impact  area  at  Fort  Greely,  Alaska  provides  an  opportunity  to  examine  the 
effects  of  munitions  use  in  an  ecosystem  as  well  as  the  ecosystem’s  response. 

The  objectives  of  this  study  are  to; 

.  Examine  the  spatial  and  temporal  characteristics  of  munitions  disturbance. 

.  Describe  and  relate  important  soil  quality  response/capacity  attributes  that 
directly  influence  the  ability  of  the  ecosystem  to  influence  the  fate  and  transport 
of  potential  UXO  contamination. 

.  Resolve  these  attributes  into  a  measurable  metric  or  indicator  that  can  then  be 
used  to  evaluate  ecosystem  response. 


Background:  Impact  Area  Design,  Disturbance  Regime  and  Ecosystem  Response 


Impact  Area  Design 

Impact  areas  are  the  ground  and  associated  airspace  within  a  training  complex 
used  to  contain  fired,  dropped  or  launched  munitions  and  any  resulting  fragments, 
debris,  and  components  (CEMML,  1999).  The  Army’s  Readiness  Program  requires  that 
impact  areas  provide  effective  individual,  crew,  and  unit  live  fire  and  maneuver 
proficiency  training  using  the  Army’s  current  and  projected  modern  weapons  systems 
(AR  21 0-21 ,  1 997,  p.  1 0).  Therefore,  impact  areas  reflect  a  mission-oriented  design 
that  is  continually  being  updated  to  accommodate  technological  and  doctrinal  changes. 

Impact  area  lands  are  more  effectively  utilized  by  grouping  firing  ranges  so  that 
safety  danger  zones  (SDZ,  i.e.,  the  hazardous,  three  dimensional  space  along  a  line 
from  the  gun  to  the  target  to  include  areas  around  the  gun  and  target)  associated  with 
individual  weapons  systems  overlap.  This  not  only  reduces  the  area  required  for  live 
fire  training,  but  also  minimizes  the  area  contaminated  by  UXO.  Additionally,  impact 
areas  are  segregated  from  surrounding  maneuver  areas  and  from  public  and  private 
lands  to  avoid  unintentional  contact  with  UXO,  which  can  be  physically  hazards.  They 
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are  intended  for  indefinite  use  (AR  210-21 ,  1997,  p.15)  because  of  the  challenges  and 
uncertainties  associated  with  clearing  UXO.  As  such,  they  represent  a  potential  source 
of  environmental  contamination  either  from  UXO  residues  or  from  other  sources  such  as 
sediment  pollution  caused  by  normal  cratering  effects.  Therefore,  current  impact  area 
design  and  management  of  older  impact  areas,  incorporates  a  high  degree  of 
environmental  planning  that  includes  pollution  prevention,  and  natural  and  cultural 
resource  protection  considerations  (TC  25-1 ,  1991 ,  p.3-4) 

Impact  area  size  primarily  depends  upon  the  size  and  type  of  units  designated  to 
use  them.  An  U.S.  Army  installation  usually  accommodates  approximately  one  Army 
division.  Divisions  typically  consist  of  1 0,000  to  1 5,000  soldiers  and  their  equipment 
and  they  have  different  missions  and  organization  such  as  light  infantry,  air  borne,  air 
assault,  mechanized  infantry  and  armor.  For  example,  armored  divisions  are  equipped 
with  hundreds  of  armored  combat  systems  such  as  tanks  and  infantry  fighting  vehicles, 
while  light  infantry  divisions  have  no  armored  combat  vehicle  systems. 

Combat  readiness  is  a  measure  of  a  unit’s  ability  to  accomplish  its  assigned 
wartime  mission.  Many  of  the  readiness  capabilities  assessed  involve  how  well 
individual  soldiers  and  combat  crews  employ  their  individual  and  crew  served  weapons 
systems.  Therefore,  Army  range  objectives  are  tailored  to  provide  ranges  and 
associated  impact  areas  to  ensure  combat  readiness  of  the  division  assigned  to  the 
installation  (Table  1 ).  Ranges  are  generally  designed  for  direct  fire  or  indirect  fire 
weapons  systems  or  for  a  combination  of  system  types.  Based  on  these  range 
requirements,  one  can  determine  an  approximate  target  and  target  array  area 
requirements  within  an  impact  area  for  a  given  set  of  ranges.  Direct  fire  weapons 
systems  such  as  tanks  and  anti  tank  systems  are  those  that  fire  to  a  point  target  such 
as  an  enemy  tank  or  fighting  position.  They  require  a  series  of  single  targets  each  of 
which  occupy  a  relatively  small  area  on  the  order  of  a  quarter  hectare.  Indirect  fire 
weapons  systems  are  those  employed  against  area  targets  such  as  a  group  of  enemy 
vehicles  or  soldiers.  They  require  groupings  of  targets  or  target  arrays  that  occupy  a 
somewhat  larger  area  on  the  order  of  one  hectare.  Typically,  division  sized  range 
facilities  and  impact  areas  require  on  the  order  of  one  or  two  thousand  targets  and 
target  arrays  occupying  many  hundreds  of  hectares  within  impact  areas. 

The  impact  area  that  is  the  subject  of  this  study  is  at  Fort  Greely  located  in 
Alaska.  Fort  Greely  has  been  an  active  Army  installation  for  over  fifty  years.  The 
installation  covers  around  267,905  hectares  and  has  34,415  hectares  of  high  hazard 
impact  areas  -  those  in  which  high  explosive  munitions  are  employed  (Figure  1  and 
Landsat  image  at  Appendix  A).  Its  impact  areas  were  designed  to  support  a  division¬ 
sized  unit,  although  today  only  one  brigade  (approximately  25%  of  a  standard  division) 
trains  at  Fort  Greely. 
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Anthropogenic  Disturbance 


Munitions  use  and  fire  resulting  from  both  munitions  use  and  maneuver  training 
are  significant  disturbance  gradients  operating  at  Fort  Greely.  An  examination  of  the 
types  of  munitions  used  in  an  impact  area,  the  quantities  expended,  and  an 
approximation  of  the  magnitude  of  potential  contaminant  loading  provides  a  measure  of 
the  problem  under  consideration.  High  explosive  munitions  are  the  focus  of  this  study 
because  they  have  both  a  destructive  and  contaminate  effect  in  the  environment.  High 
explosive  munitions  are  those  that  explode  on  or  near  the  ground  and  therefore  have 
the  potential  of  creating  a  crater.  Other  munitions  include  smoke  projectiles,  illumination 
projectiles,  and  small  arms  and  inert  projectiles  (CEMML,  1999).  These  munitions  have 
the  potential  of  creating  craters,  but  this  potential  is  relatively  small  compared  to  high 
explosive  munitions. 

The  magnitude  of  the  physical  disturbance  occurring  within  Fort  Greely  impact 
areas  can  be  approximated  by  an  evaluation  of  munition  expenditure  data  for  the  period 
1992-1997  (Table  2  and  Tabular  Data  for  Munitions  Loading  at  Appendix  B).  Fort 
Greely’s  high  hazard  impact  areas  had  at  least  185,846  high  explosive  munitions  fired, 
launched  or  dropped  into  its  impact  areas  or  approximately  30,974  munitions  per  year 
(CEMML,  1999).  To  account  for  the  area  exposed  to  high  explosive  munitions,  one  can 
approximate  the  surface  area  of  target  arrays.  This  area  could  encompass  up  to  945  ha 
exposed  to  munitions  disturbance,  or  about  2.75%  of  the  total  impact  area  lands. 
Determining  the  approximate  annual  cratered  surface  area  refines  this  estimate  of 
disturbance.  Again,  using  the  CEMML  (1 999)  data  set  and  assuming  a  crater  radius  for 
each  type  of  munition  it  can  be  determined  that  approximately  10.21  ha  were  disturbed 
annually  during  this  period  representing  approximately  1 .08%  of  the  lands  exposed  to 
munitions  or  less  than  0.03%  of  the  total  high  hazard  impact  areas. 

In  addition  to  the  inherent  kinetic  and  blast  destructive  effects  of  munitions,  there 
is  a  little  understood  contaminant  effect  associated  with  munitions  that  fail  to  function 
properly.  There  are  two  modes  of  munition  failure  —  dud  and  low  order  detonation.  A 
dud  is  a  round  that  completely  fails  to  function  at  the  target  while  a  low  order  detonation 
is  a  round  that  only  partially  functions  leaving  part  of  the  high  explosive  filler  scattered 
over  the  target  area  (USATCES,  2000).  Dud  rounds  retain  their  explosive  filler  inside 
the  delivery  body  of  the  munition  and  are  relatively  isolated  from  the  environment. 
However,  during  impact  the  delivery  body  can  become  fractured  and  over  time  there  is 
the  possibility  that  small  amounts  of  explosive  filler  can  be  mobilized  through  dissolution 
with  water.  Low  order  detonations  can  result  in  significant  amounts  of  ready  availability 
explosive  filler  for  dissolution  and  subsequent  movement  into  the  environment  (Brannon 
et  al.,  1999).  Therefore,  both  dud  and  low  order  detonation  munitions  have  the  potential 
for  environmental  contamination  by  leached  explosive  filler. 

Explosive  fillers  such  as  Trinitrotoluene  (TNT),  Royal  Demolition  Explosive 
(RDX),  and  High  Melting  Explosive  (HMX)  are  subject  to  environmental  fate  and 
transport  congruent  with  the  ecological  setting  in  to  which  they  are  released.  Release  is 
contingent  on  dissolution  rates  and  there  have  been  a  limited  but  varied  number  of 
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reported  rates.  Brannon  et  al.  (1999)  reported  the  dissolution  fluxes  for  free  solid  to 
liquid  phase.  Their  model  considers  dissolution  flux  critical  for  quantifying  contaminant¬ 
loading  rates.  They  report  TNT  has  the  highest  rate  at  4,164  pg  cm  '  hr  HMX  has  a 
rate  of  454  pg  cm  hr  and  RDX’s  has  the  lowest  dissolution  rate  at  360  pg  cm  hr 
Lynch  et.  al,  (2001)  report  much  higher  rates  of  604,800  pg  cm  hr  25,056  pg  cm 
hr  and  8294  pg  cm  hr  for  TNT,  HMX  and  RDX  respectively.  They  used  weapons 
grade  explosives  and  subjected  samples  to  mixing,  which  may  be  more  similar  to 
particular  environmental  conditions  such  as  free  product  exposed  to  surface  water  flow. 
Since  dissolution  is  a  function  of  interfacial  surface  area,  temperature,  and  energy  input 
(i.e.  degree  of  fluid  motion),  they  may  have  observed  these  higher  rates  due  to  relatively 
large  energy  input.  Also,  dissolution  rates  of  pure  explosive  will  differ  from  the  rate  of 
mixtures  of  explosives  used  in  munitions  (e.g.,  an  artillery  munition  may  have  a  mixture 
of  TNT  and  RDX).  This  is  because  most  manufacturing  processes  result  in  one 
explosive  encapsulated  within  another  explosive,  or  within  compounds  that  serve  as 
explosive  binders.  For  example,  “certain  formulations  include  binders  made  of  insoluble 
wax;  explosives  almost  entirely  encapsulated  in  wax  would  persist  for  extended  periods 
of  time  yet  yield  comparably  lower  concentrations.”  (Lynch,  2001 ,  personal 
communication). 

The  magnitude  of  the  contaminant  threat  can  be  approximated  by  an  evaluation 
of  munition  failure  rates.  Using  failure  rates  published  by  USATCES  (2000)  and 
assuming  the  same  failure  rates  are  applicable  to  Air  Force  munitions,  analysis 
indicates  that  there  could  be  approximately  1 ,271  dud  and  95  low  order  munitions 
failures  per  year  on  Fort  Greely  impact  areas.  Loading  rates  could  be  derived  from  this 
data,  but  the  highly  variable  nature  of  the  state  of  explosive  fillers  available  for 
mobilization  (e.g.,  fully  contained  in  a  dud  delivery  body,  partially  contained  in  a  low 
order  delivery  body  and/or  scattered  by  a  low  order  detonation)  coupled  with  varying 
solubility  and  dissolution  fluxes  resulting  from  manufacturing  processes  makes  any 
attempt  very  uncertain.  The  most  conservative  approach  to  assume  that  all  of  the 
explosive  filler  from  dud  and  low  order  munitions  failures  is  available  for  dissolution. 
Assuming,  optimum  dissolution  conditions  this  results  in  a  maximum  potential  annual 
average  of  4,460  kg  of  high  explosive  materials  and  translates  into  a  maximum 
theoretical  loading  rate  within  the  target  array  area  exposed  to  high  explosive  munitions 
of  about  0.46  gm  m'^  yr  Alternatively,  assuming  that  this  maximum  potential  annual 
average  of  high  explosive  materials  is  deposited  in  cratered  areas,  the  maximum 
theoretical  loading  rate  could  be  as  high  as  44  gm  m'^  yr  (Table  2).  Assuming  this 
load  resides  in  the  top  1  cm  of  the  soil  profile,  a  soil  bulk  density  of  2.0  Mg  m'^  (e.g., 
glacial  till;  Brady  and  Weil,  1999,  p.  136),  and  that  no  transformations  take  place,  this 
equates  to  potential  concentrations  ranging  from  0.046  pg  gm'^  for  target  areas  to  4.4 
pg  gm'^  for  craters  for  a  given  year.  This  is  potentially  well  below  Environmental 
Protection  Agency  residential  risk  based  concentrations  (i.e.,  5.8  to  3900  pg  gm'^  for 
various  explosives)  reported  by  Walsh  et  al.,  (2001). 

Finally,  the  other  major  form  of  human  disturbance  is  in  the  form  of  fires  started 
either  by  maneuver  training  or  by  munitions  use.  Jorgenson  et  al.  (2001)  report  that  fire 
frequencies  and  percent  area  burned  are  slightly  higher  than  similar  forest  types  in 
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interior  Alaska.  Since  1950, 1.2%  of  Fort  Greely  has  burned  annually  whereas  interior 
Alaska  forestland  has  exhibited  percent  annual  burn  areas  ranging  from  1%  between 
1940  and  1970,  to  6  %  between  1970  and  the  present.  This  slightly  higher  percentage 
is  consistent  with  the  effects  of  military  training  (Demarais  et  al.,  1999).  Fire  is 
significant  in  this  landscape  because  it  is  the  dominant  process  affecting  succession 
(Jorgenson  et  al.,  2001)  and  is  ultimately  important  in  any  soil  forming  process. 


Soil,  Climate  and  Hydrologic  Ecosystem  Component  Response 

When  considering  ecosystem  response  in  the  face  of  disturbance  one  naturally 
thinks  of  the  manner  in  which  a  disturbed  system  recovers.  A  measure  of  stability  is  the 
concept  of  resilience  or  the  time  it  takes  a  system  to  return  to  a  steady  state  following 
perturbation  (DeAngelis  et  al.,  1989).  Climate,  soils,  hydrology,  flora  and  fauna  interact 
and  affect  energy,  moisture,  nutrient  and  disturbance  gradients  which  in  turn  regulate 
the  structure  and  function  of  ecosystems  and  their  ability  to  rebound  from  disturbances 
(Jorgenson  et  al,  2001).  High  system  resilience  has  been  characterized  in  terms  of  high 
energy  fluxes  through  food  webs  and  short  nutrient  cycling  times  (O’Neill,  1976; 
DeAngelis,  1980,  1989;  Harwell  et  al.,  1977;  and  Loreau,  1994).  In  the  case  of  impact 
areas,  expected  soil  component  responses  might  be  recovery  from  cratering  through 
organic  matter  inputs,  and  the  assimilation  and  transformation  of  explosive  compound 
contamination  within  the  soil  matrix.  These  responses  should  be  intimately  tied  to 
biomass  inputs,  and  to  processes  that  influence  nutrient  availability  within  the  soil 
component. 

Fate  and  transport  processes  are  continuously  at  work  within  the  environment. 
Brannon  et  al.  (1999)  developed  a  conceptual  model  for  fate  and  transport  of  UXO  that 
provides  a  basis  for  a  qualitative  evaluation  of  the  key  processes  at  work  in  the 
environment.  Their  model  outlined  representations  for  soil,  sediment,  and  water 
environments,  described  source  term  considerations,  and  developed  descriptor 
formulations  of  key  processes  resulting  in  explosive  contaminant  mobilization, 
retardation,  and  transformation  effects.  With  modifications,  these  key  processes  are 
dissolution,  the  initiating  release  mechanism  described  above;  adsorption  and  chelation 
processes  in  soils,  considered  to  be  an  important  retardation  mechanism;  and  biotic 
transformation,  or  for  the  purposes  of  this  study  biological  action,  consisting  of 
metabolism  by  microorganisms  and  uptake  by  plants.  Biological  action  is  regarded  as 
an  important  transformation  mechanism  (Comfort  et  al.,  1995;  Larson  et  al.,  1998a; 
Pennington,  1988;  Price  et  al.,  1998;  Brannon  et  al  1999;  Pennington  et  al.,  1999a;  and 
Pennington  et  al.,  1999b). 

Applying  this  conceptual  model  into  an  Alaskan  ecological  framework  involves 
relating  process  and  effect  to  relevant  ecosystem  components  and  their  governing 
ecological  setting.  For  example,  dissolution  is  the  process  by  which  TNT,  RDX  and 
HMX  are  dissolved  by  a  solvent  from  free  product  to  solutes.  In  this  case,  the  solvent  is 
water  that  is  provided  via  the  hydrologic  cycle.  The  hydrologic  cycle  is  an  ecological 
component  that  determines  the  availability  of  water,  and  its  relative  importance  is 
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governed  by  the  ecological  setting  of  the  region  in  which  UXO  is  located.  In  Fort 
Greely,  water  is  available  in  an  unfrozen  form  for  relatively  short  time  frames. 

Adsorption  is  the  process  by  which  solutes  move  to  and  bind  with  a  sorbent. 

This  processes  is  reversible  and  is  dependent  on  the  partitioning  affinity  of  the  solute  for 
either  the  solute  or  sorbent  (LaGrega  et  al.,  1994).  For  this  analysis,  the  process  is  the 
movement  or  partitioning  of  dissolved  TNT,  RDX,  and  HMX  to  soil  colloidal  surfaces. 
Generally,  soil  colloidal  surfaces,  composed  of  clays  and  humic  matter,  are  negatively 
charged.  In  the  subsurface,  adsorption  is  primarily  a  process  of  positively  charged 
organic  chemicals  accumulating  at  these  colloidal  soil  surfaces  of  which  naturally 
occurring  humic  matter  in  soil  plays  a  significant  role  (LaGrega  et  al.,  1994).  Therefore, 
the  notable  parameter  for  adsorption  is  colloidal  content  of  the  soil  component.  Colloid 
content  can  be  measured  via  a  soil’s  cation  exchange  capacity  (CEC)  which  in  turn  is 
used  to  define  a  partitioning  coefficient  for  contaminant  to  soil  adsorption.  In  terms  of 
the  influence  of  ecological  setting  on  the  adsorption  process,  one  can  evaluate  the  soil 
quality  and  its  CEC  to  obtain  an  indication  of  the  relative  retardation  potential  of  the  area 
in  which  UXO  is  found.  Generally,  Alaskan  soils  are  not  as  weathered  as  soils  in  more 
temperate  regions.  Therefore  they  generally  have  less  clay  fraction  and  humus 
becomes  the  dominant  contributor  to  CEC. 

Adsorption  of  TNT  to  soil  has  been  documented  and  considered  to  be  very 
important  in  immobilizing  TNT  and  its  transformation  products  (Larson  et  al.,  1998b; 
Pennington  et  al.,  1999c:  Comfort  et  al,  1995;  Weissmahr  et  al.,  1998;  Myers  et  al., 

1998)  and  in  particular  to  humic  soils  (Hundal  et  al.,  1997).  TNT  adsorption  has  been 
positively  correlated  with  cation  exchange  capacity  (CEC)  of  soils  (Pennington,  1988; 
Brannon  et  al.,  1999).  RDX  has  been  shown  to  sorb  at  lower  rates  than  TNT  (Singh  et 
al.,  1998;  Price  et  al.,  1998)  and  also  can  be  positively  correlated  with  soil  CEC 
characteristics  (Brannon  et  al.,  1999).  While  HMX  adsorption  occurs,  little  is  known 
about  the  solute-sorbent  process  to  enable  the  prediction  of  adsorption  coefficients  for 
HMX  (Brannon  et  al.,  1999). 

Metals  availability  are  largely  determined  by  the  colloidal  portion  of  soil  organic 
matter  (Tate,  1987).  Humic  acids  contain  highly  reactive  anionic  groups  that  chelate  or 
complex  metals.  Tate  (1987)  states  that  due  to  the  high  reactivity  of  these  groups  there 
is  a  positive  relationship  between  degree  of  humification  and  complexation.  This 
relationship  should  be  less  pronounced  in  acidic  soils  as  may  be  the  case  in  Alaska 
were  coniferous  forests  intensify  soil  acidity. 

Biological  action  consists  of  the  mechanisms  of  biodegradation  and  bio-uptake  of 
contaminants.  Biodegradation  involves  the  oxidation  and  reduction  of  organic 
contaminants  and  bio-uptake  of  orgainic  and  inorganic  contaminants  involves  plant  root 
uptake  by  vegetation.  Both  depend  on  the  ecological  components  of  soil  and  climate. 
Soil  structure  provides  habitat  for  microorganisms,  which  obtain  energy  and  carbon  from 
the  soil.  Vegetation  requires  soil  for  structural  support,  for  a  source  of  oxygen  for  root 
respiration,  and  for  a  source  of  macro  and  micronutrients.  Soil  microorganisms  and 
plants  both  require  water  and  a  temperature  regime  conducive  for  growth,  which  is  a 
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function  of  climate.  Soil  organic  carbon  (SOC)  content  and  climate  in  terms  of  effective 
growing  season  (EGS)  are  good  indicators  of  conditions  favorable  for  microbial  and 
plant  activity.  SOC  and  EGS  can  provide  a  relative  measure  of  the  transformation 
potential  of  the  region  in  which  UXO  is  found.  In  Alaska,  growing  seasons  are  relatively 
short  so  biological  activity  may  not  be  as  great  as  in  more  temperate  regions. 

Biologically,  TNT  can  undergo  both  aerobic  and  anaerobic  degradation  (Brannon 
et  al  1999;  Larson  et  al.,  1998a;  Pennington  et  al.,  1999a)  provided  there  are  sufficient 
carbon  and  energy  sources  for  microorganisms  (Pennington  et  al.,  1 999b).  RDX  and 
HMX  undergo  biodegradation  only  under  anaerobic  conditions  (Price  et  al.,  1998; 
Comfort  et  al.,  1995  as  reported  by  Larson  et  al.,  1998a).  Vegetation  also  can  act  as  a 
sink  for  explosive  contaminants  by  root  uptake  (Larson  et  al.,  1998a). 

All  of  the  explosives,  as  with  any  solute,  can  be  transported  by  groundwater 
advection.  These  UXO  fate  and  transport  characterizations  help  to  serve  as  the 
conceptual  foundation  for  the  development  of  a  framework  for  the  ecological 
assessment  of  impact  area  resilience. 


Spatial  and  Temporal  Scaling  of  Ecological  Soil  Evaluation 

An  ecologicai  soil  evaluation  is  an  examination  of  soil  quality  and  of  ecosystem 
response  to  disturbance.  Therefore,  this  research  effort  is  designed  based  on 
appropriate  scaling  of  the  soil  component  processes  and  of  disturbance  response. 

Also,  as  one  begins  to  evaluate  the  ability  or  capacity  of  an  impact  area  ecosystem  to 
respond  to  munitions  disturbance  one  must  be  concerned  with  the  relevant  spatial  and 
temporal  scales  associated  with  impact  area  use.  Therefore,  the  first  step  in  this 
evaluation  process  is  one  that  determines  the  appropriate  spatial  and  temporal  scales  of 
interest  relative  to  munitions  disturbance  and  UXO  contamination. 

O’Neill  et  al.  (1986)  argue  that  a  critical  step  in  addressing  any  problem  is 
determining  the  scale  of  observation.  This  allows  one  to  focus  on  a  particular 
organizational  level  characterized  by  a  range  of  rates  that  can  be  segregated  from 
higher  and  lower  levels  of  dynamics  thereby  describing  a  system  that  is  conducive  to 
study  the  results  of  which  are  not  confounded  by  indifferent  dynamics.  This  approach 
is  further  supported  by  Wiens  (1989)  who  states  that  the  goal  of  an  investigation  should 
be  to  obtain  a  balance  between  extent  and  grain  (i.e.,  the  upper  and  lower  limits  of 
resolution  respectively)  in  order  to  minimize  the  spatial  variance  of  a  process  of  interest. 
Ecological  setting,  munitions  use  effects,  soil  quality  indicators  and  information 
management  requirements  are  scalar-sieves  that  can  be  used  to  define  appropriate 
scales  of  interest  and  relevant  parameters.  Ecological  setting  provides  a  general 
indication  for  appropriate  upper  level  scaling  rational.  An  examination  of  munitions  use 
within  an  impact  area  helps  to  define  the  anthropogenic  disturbance  regime  resulting  in 
a  higher  degree  of  spatial  and  temporal  resolution.  Soil  quality  indicators  help  to 
describe  relevant  large-scale  (i.e.,  local)  processes  and  provide  a  refined  indication  for 
both  spatial  and  temporal  scaling  in  terms  of  ecosystem  response.  Finally,  information 
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management  requirements  integrate  all  of  these  scalars  into  an  investigative 
arrangement. 


Ecological  Settings 

The  ecological  settings,  or  the  dominant  ecosystem  components  of  interest,  are 
important  in  helping  to  determine  spatial  and  temporal  scales  of  interest.  As  described 
above,  climate,  soils,  hydrology,  flora  and  fauna  are  examples  of  ecosystem 
components.  They  interact  and  they  affect  energy,  moisture,  nutrient  and  disturbance 
gradients  which  in  turn  regulate  the  structure  and  function  of  ecosystems  (Jorgenson  et 
al,  2000)  found  within  the  boundaries  of  an  impact  area.  They  involve  complex 
processes  that  demand  a  hierarchical  approach  designed  to  define  the  system  of 
interest  and  the  types  of  measurements  and  their  interpretations  taken  on  the  system 
(O’Neill  et  al.,  1986).  This  allows  one  to  determine  and  focus  on  a  component  of 
interest  and  its  range  of  variation.  Components  that  have  slower  rates  become 
background  and  components  with  faster  rates  are  averaged  out  (O’Neill  et  al.,  1986). 
Relevant  components  in  this  evaluation  are  climate,  flora  and  soils. 

The  climate  of  the  system  determines  the  availability  of  liquid  water,  which  will 
largely  determine  the  potential  for,  contaminate  mobilization  (Brannon  et  al.,  1999). 
Also,  it  is  a  major  determinant  in  vegetative  growth  and  soils  development.  However,  it 
operates  at  spatial  temporal  scales  that  are  at  higher  levels  and  at  slower  rates  of 
change  than  the  disturbance  response  phenomena  of  interest.  Climate  is  therefore 
considered  a  background  component. 

Since  the  majority  of  the  contaminant  threat  lies  at  or  near  the  surface  of  the 
ground,  the  soil  component  is  an  important  ecological  setting  for  an  impact  area  at  Fort 
Greely.  The  flora  component  and  the  successional  characteristics  of  the  system  are 
intimately  tied  to  soil  development  and  therefore  soil  component  processes  should  be 
reflected  at  the  forest  stand  level.  If  one  considers  the  extent  of  the  study  to  be  the 
boundaries  of  the  impact  area  and  the  major  forest  ecotypes  that  populate  It,  then  the 
grain  should  be  at  a  scale  that  minimizes  spatial  variance  within  these  boundaries.  With 
extent  held  as  a  constant  scale  would  be  held  constant  the  grain  would  be  scaled  at  the 
level  of  sub  forest  stand.  This  should  be  large  enough  to  decrease  between  grain 
spatial  variance  within  the  extent  of  the  study  area. 


Munitions  Use  in  Impact  areas  -  Levels  of  Disturbance 

Walker  and  Walker  (1991),  describe  a  hierarchy  of  disturbance  scales  in  Arctic 
Alaska.  They  state  that  anthropogenic  disturbances  are  generally  on  the  microscale 
level  (10'^  to  10®  m^),  and  that  cumulative  impacts  are  on  the  mesoscale  level  (10®  to 
10^®  m^).  This  is  particularly  true  of  munitions  use  in  impact  areas  wherein  munitions 
cause  craters  that  range  in  size  from  less  than  a  meter  to  several  meters  in  diameter. 
Also,  impact  areas  accommodate  firing  of  weapon  systems  that  have  ranges  on  the 
order  of  tens  of  kilometers  and  must  retain  the  effects  munitions  requiring  several 
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hundred  to  thousands  of  hectares.  Temporally,  munitions  use  within  impact  areas  falls 
within  the  natural  disturbance  frequency  range  between  1  and  1000  years  (e.g.,  monthly 
use  over  a  period  of  fifty  years),  and  is  within  the  same  recovery  timeframe.  Therefore, 
an  evaluation  that  is  intended  to  examine  munitions  use  over  spatial  and  temporal 
scales  of  hectares  to  square  kilometers  and  years  to  hundreds  of  years  respectively, 
falls  comfortably  into  a  microscale  to  mesoscale  spatial  and  temporal  realm  as  defined 
by  Walker  and  Walker  (1991).  Finally,  these  scales  of  observation  lend  themselves  to 
data  derived  from  satellite  sensed  data  and  surface  surveys  to  examine  anthropogenic 
effects  such  as  cratering  and  ecosystem  recovery. 


Inherent  Soil  Quality  Indicators 

Soil  quality  should  be  affected  by  munitions  use  primarily  through  physical 
disturbance.  Contaminant  fate  and  transport  in  soils  is  dependent  on  same  soil 
chemical  properties  that  influence  nutrient  availability  and  uptake  by  plants  and 
microorganisms  (Heil  and  Sposito,  1997).  Therefore,  an  examination  of  significant 
physical,  chemical  and  biological  indicators  and  their  associated  processes  will  provide 
an  indication  of  appropriate  spatial  and  temporal  scales  of  study. 

Topp  et  al.,  1997  describe  several  physical  attributes  of  soil  quality.  Key  among 
them  is  soil  water  transmission  and  soil  strength.  Soil  water  transmission  capability 
determines  the  rate  at  which  excess  soil  water  drains  out  of  the  root  zone  and  reflects 
system  water  availability  and  matter  transfer  potentials.  It  is  described  by  the  soil’s 
saturate  hydraulic  conductivity  with  rates  ranging  from  1 0^  to  1 meters/year  for  fine  to 
coarse  textured  soils  respectively  (Topp  et  al.,  1997).  Soil  strength  is  a  function  of 
cohesion  and  friction  and  is  a  soil’s  ability  to  resist  shear  breakage,  which  is  potentially 
important  in  physical  disturbance  response.  It  can  be  assessed  over  short  spatial  and 
temporal  scales  of  square  meters  and  years.  Both  soil  water  transmission  and  strength 
can  be  related  to  soil  texture  (Topp  et  al.,  1997). 

Chemical  attributes  of  importance  include  mineralogy,  pH,  buffer  capacity, 
organic  matter,  and  cation  exchange  capacity.  They  affect  rates  of  reaction  for 
processes  such  as  soluble  complex  formation,  adsorption,  mineral  dissolution,  are 
relevant  to  plant  growth,  (i.e.  available  nutrient  utilization,  mobility  of  nutrient  species, 
and  diffusion  processes)  and  are  on  the  order  of  meters/year  (Heil  and  Sposito,  1997). 

A  notable  biological  input  into  the  chemical  system  is  humus  in  organic  matter  provides 
reactive  colloidal  adsorption  sites  for  nutrient  storage  and  in  doing  so  affects  cation 
exchange  capacity  (Heil  and  Sposito,  1997).  Organic  matter  has  transformation  and 
colloidal  production  rates  that  range  from  a  few  years  to  centuries  (Parton,  1996). 

Gregorich  et  al.  (1997)  describe  two  primary  biological  indicators  of  soil  quality. 
They  are  total  organic  carbon  and  nitrogen  and  are  related  to  soil  organic  matter 
dynamics  resulting  from  management  practices  such  as  tillage.  Cratering  effects  of 
munitions  can  be  related  to  tillage.  Long  term  studies  are  required  to  observe 
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significant  changes  soil  organic  matter  accumulation  or  loss.  Associated  spatial  and 
temporal  scales  are  on  the  order  of  hectares  and  tens  to  hundreds  of  years. 

Overall,  spatial  and  temporal  scales  of  interest  are  bounded  by  physical, 
chemical  and  biological  properties  of  the  soil  component.  Chemical  properties  seem  to 
be  generally  aligned  with  small  spatial  and  temporal  time  scales  relevant  to  plant 
growth.  Physical  properties  are  generally  aligned  with  slightly  larger  scales.  Biological 
indicators  are  somewhat  larger  than  chemical  and  physical  indicators,  and  on  the  order 
of  hectares  and  tens  to  hundreds  of  years.  Also,  biological  indicators  are  significant  in 
that  they  represent  a  higher  level  of  system  dynamics  controlling  chemical  and  physical 
soil  properties.  Based  on  the  discussion  above,  key  indices  for  soil  quality  as  they 
relate  to  munitions  use  may  be: 

.  pH  and  Buffer  Capacity 
.  Cation  Exchange  Capacity 
.  Nutrient  Status  (e.g.,  N,  P,  K) 

.  %  Sand,  Silt,  Clay  (Texture) 

.  %  Organic  Carbon  and  Organic  Matter 


Information  Management 

Information  management  allows  one  to  understand  technological  limitations  and 
capabilities  and  at  the  same  time  helps  to  integrate  relevant  spatial  and  temporal  scales 
into  a  system  which  is  not  confounded  by  higher  and  lower  levels  of  dynamics.  In 
effect,  it  allows  us  to  determine  sampling  methodology  and  subsequent  analysis. 

Sample  plot  size  is  a  compromise  between  the  area  utilized  by  a  typical  target  array,  the 
spatial  resolution  of  sensor  and  mapping  software,  and  an  area  that  is  conducive  to 
effective  land  management.  A  one-hectare  plot  size  meets  all  of  these  conditions. 
Target  arrays  range  in  size  from  100  m^  for  direct  fire  targets  to  1  -4  ha  for  indirect  fire 
targets:  one  hectare  equates  to  12.3  pixels  of  Landsat  TM  data  (28.5  m  spatial 
resolution);  filtered  (i.e.,  a  3x3  or  a  5x5  majority  filter)  thematic  classifications  of  TM  data 
resolve  to  approximately  0.7-2  ha;  and  one  hectare  is  a  universally  used  management 
unit.  (Figure  2).  Temporally,  remotely  sensed  data  can  provide  recent  and  timely 
information  about  disturbance  and  response.  Finally,  soil  quality  information  can  be 
easily  obtained  and  readily  assessed  at  these  scales. 

Ecological  settings  provide  a  macro  scale  perspective,  which  is  critical  to 
understanding  higher  lever  background  dynamics.  Munitions  use  helps  to  further  define 
the  extent  of  the  problem  narrowing  the  range  of  spatial  and  temporal  scales.  Soil 
quality  is  dominated  by  process  dynamics  that  operate  at  scales  relevant  to  the  grain  of 
measure  and  represents  a  micro  perspective.  Information  management  allows  one  to 
realize  a  balance  between  higher  and  lower  level  of  dynamics  within  technological 
capability. 
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Methods 


Study  and  Control  Area  Descriptions 

Fort  Greely  Military  Installation  is  located  in  central  Alaska  approximately  160  km 
to  the  southeast  of  Fairbanks,  near  the  city  or  Delta  Junction.  The  installation  covers 
approximately  267,900  ha  and  consists  of  the  cantonment  area.  West  Training  Area  and 
East  Training  Area.  Fort  Greely  is  a  US  Army  training  and  testing  installation  with 
approximately  34,400  ha  of  high  hazard  impact  areas.  The  Washington  impact  area  is 
the  subject  of  this  study  and  is  used  for  a  variety  of  weapons  testing  and  training 
involving  mortar  and  artillery  firing  exercises,  tank  gunnery,  and  aerial  bombing. 
Washington  impact  area  is  set  in  a  rivierine  floodplain  and  encompasses  approximately 
800  ha  within  the  Delta  River  floodplain  (Walsh,  et  al.,  2001).  At  the  south  boundary  of 
the  impact  area  is  a  portion  of  the  Washington  Range  (i.e.,  an  area  from  which  weapons 
systems  are  fired  into  the  impact  area).  This  area  extends  approximately  2.5  kms  to  the 
south,  is  very  similar  to  the  impact  area  and  is  a  control  area  for  the  study.  It  is  called 
the  Washington  Range  control  area  (Figure  1). 

The  second  control  area  Jarvis  Creek  control  area  is  approximately  15  km  to  the 
east-north-east  of  the  Washington  impact  area  is  used  for  maneuver  training  and  is 
named  for  Jarvis  Creek  that  runs  along  the  training  area’s  western  boundary.  The  area 
is  similar  to  Delta  River  in  terms  of  direction  of  water  flow,  elevation,  geomorphic  units 
and  ecotypes  (Jorgenson  et  al.,  2001). 

Further  description  of  control  areas  and  study  area  can  be  made  in  terms  of  the 
dominant  ecosystem  components  of  the  Fort  Greely  area,  which  are  climate, 
physiography,  hydrology  and  disturbance.  The  climate  of  Fort  Greely  as  described  by 
Bailey  (1998)  is  Boreal  Subarctic  characterized  by  a  great  seasonal  range  in 
temperatures.  While  its  climate  is  moist  all  year.  Fort  Greely  has  relatively  small  annual 
amounts  of  precipitation  that  is  concentrated  in  warm  months.  Winter  is  the  dominant 
season  with  cool  short  summers.  Weather  data  covering  63  years  from  Big  Delta 
FAA/AMOS  AP,  AK  (500770)  indicates  that  January  is  the  coldest  month  with  an 
average  maximum  temperature  of  -16.1C  and  average  minimum  temperature  of  - 
24. 1C.  July  is  the  warmest  month  with  an  average  maximum  temperature  of  20.8C  and 
average  minimum  temperature  10.3C.  Total  average  annual  precipitation  is  29.5cm 
with  July  having  the  highest  monthly  average  of  6.7cm.  Total  average  annual  snowfall 
amounts  to  1 1 1 .3  cm  with  October  having  highest  monthly  average  snowfall  of  23.4cm. 
Average  snow  depth  during  winter  is  10.2  cm  with  the  month  of  February  having  the 
highest  monthly  average  depth  of  25.4  cm.  Using  the  USACERL  methodology  (Diersing 
et  al.,  1 990),  Fort  Greely’s  effective  growing  season  is  25  May  -  31  August  (Figure  3). 

The  region’s  physiography  is  described  as  the  Tanana-Kuskokwim  Lowlands 
(CEMML,  1999).  It  consists  of  a  broad  depression  north  of  the  Alaska  Range  on  top  of 
coalescing  outwash  fans  that  spread  northward.  The  area  is  generally  underlain  by 
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permafrost,  however  areas  over  deep  gravels  and/or  on  south  facing  slopes  are 
relatively  free  of  permafrost.  Alluvium  consists  of  porous  gravel  sands.  The  dominant 
geomorphic  feature  associated  with  the  impact  area  and  control  areas  is  glaciofluvial 
sediment  deposits  from  the  past  Delta  Glaciation  and  from  modern  glaciers  (CEMML, 
1999).  Flooding  events  are  generally  infrequent,  as  these  areas  are  located  primarily 
on  inactive  glacial-fluvial  outwash  deposits  or  in  inactive  floodplains  of  meandering  and 
headwater  streams.  Soils  are  poorly  developed  Inceptisols,  underdeveloped  Entisols 
and  Histosols.  Both  stream  systems  physiographies  are  classified  as  riverine  and  flat 
(Jorgenson  et  al.,  2001). 

The  hydrology  of  both  areas  is  similar  and  characterized  by  braided  glacial 
streams  over  glacial  outwash  fans.  Both  the  Delta  River  and  Jarvis  Creek  flow 
generally  northward  from  Alaska  Range  to  Tanana  River.  Delta  River,  which  runs 
through  Washington  Impact  Area  and  Fort  Greely,  drains  4305  square  kilometers  and  is 
130  kilometers  long.  It  primary  source  is  glacial  melt  water  from  the  Alaska  Range  and 
is  entrenched  up  to  60  meters  into  an  alluvial  fan.  The  Jarvis  Creek,  which  runs  through 
control  area  and  Fort  Greely,  drains  96  square  kilometers,  is  65  kilometers  long  and  is 
entrenched  into  same  alluvial  fan.  It  is  an  influent  stream  that  is  dry  during  the  winter, 
and  has  high  water  during  spring  melt  and  summer  rains  (CEMML,  1999).  Delta  River 
and  Jarvis  Creek  are  perched  above  single  aquifer  with  varying  local  confinement. 

They  both  streams  lose  some  flow  to  the  aquifer  in  winter.  The  groundwater  table  is 
less  than  10  meters  below  Delta  River  and  flows  to  the  northeast  (CEMML,  1999). 
Finally,  the  natural  disturbance  component  for  both  areas  consist  primarily  of  fluvial  and 
fire  processes  and  wind  erosion  (Jorgenson  et  al.,  2001). 


Ecotype  Descriptions 

Jorgenson  et  al.,  (2001)  is  the  primary  source  for  ecological  description  of  Fort 
Greely.  They  utilized  a  hierarchical  classification  scheme  to  present  three  levels  of 
ecosystem  organization  and  descriptions  for  the  Fort  Greely  installation,  which  are: 

.  Ecodistrict  -  (1 :500,000  scale)  Middle  Tanana  Floodplain  Common  Regional 
Ecosystem  with  the  Delta  River  Floodplain,  Jarvis  Creek  Floodplain  sub  districts. 

.  Ecosection  -  (1 :1 00,000  scale)  Glacial  Fluvial  Outwash  -  Active-riverbed  Deposit, 
Inactive-riverbed  Deposit,  Abandoned  Riverbed,  and  Inactive  Cover  Deposit. 

.  Ecotype-  (1 :25,000  -  50,000  scale)  areas  within  Fort  Greely  with  homogeneous 
topography,  terrain,  soil,  surface  form,  hydrology  and  vegetation. 

Their  ecotype  descriptions  are  is  based  primarily  on  physiography  and  vegetation 
structures  with  the  vegetation  type  or  successional  stage  the  differentiating 
characteristics.  They  argue  that  this  approach  is  effective  at  differentiating  soil 
characteristics  and  vegetation  composition  and  that  the  integration  of  these 
characteristics  is  useful  in  discriminating  ecotypes  that  “may  have  different  sensitivities 
to  disturbance.”  (Jorgenson  et  al.,  2001 ,  p.18)  which  is  of  particular  interest  for  this 
evaluation. 
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The  physiography  for  both  areas  described  above  are  classified  as  riverine. 

There  are  five  riverine  ecotypes  associated  with  the  study  and  control  areas.  These  are 
differentiated  by  their  vegetative  cover  and  consist  of  scrub,  broadleaf  forest,  mixed 
forest  and  needleleaf  forest  ecotypes  and  are  described  below. 

The  scrub  ecotype  is  what  Jorgenson  et  al.,  (2001)  describe  as  Riverine  Dwarf,  Low 
and  Tall  Scrub.  Vegetation  consists  of  a  range  of  early  successional  scrub  species  that 
are  highly  variable  and  frequently  includes:  Dryas  drummondii,  Populus  balsamifera 
saplings,  Shepherdia  canadensis,  Fragaria  virginiana,  Oxytopis  campestris,  Elaeagnus 
commutata,  Potentilla  multifida,  Salix  alaxensis,  Potentella  fruticosa,  Alnus  tenuifolia, 
Salix  planifolia,  Betula  nana,  and  Calamagrostis  canadensis.  Soil  profiles  range  from 
stratified  to  massive  gravel  with  occasional  thin  and  sandy  layers  that  lack  organics  to 
interbedded  silts  and  sands  with  thin  surface  organic  layers.  These  soils  are  generally 
well  drained  (i.e.,  moist)  to  excessively  drained  (i.e.,  dry)  and  exhibit  activities  that  range 
from  slightly  alkaline  to  slightly  acidic. 

The  broadleaf  forest  ecotype  is  what  Jorgenson  et  al.,  (2001)  describe  as  Riverine 
Broadleaf  Forest.  Vegetation  consists  of  a  range  of  early  to  mid  successional  broadleaf 
species.  Open  to  closed  canopies  are  composed  of  Populus  balsamifera  (occasionally 
mixed  with  Populas  tremuloides),  and  Picea  glauca.  Understories  consists  of 
Shepherdia  canadensis,  Potentella  fruticosa,  Dryas  drummondii,  Fragaria  virginiana, 
Ceratodon  purpureus.  Astragalus  spp.,  Geocaulon  lividum  and  Linnaea  borealis.  Soil 
profiles  range  from  interbedded  gravel,  sand  and  silt,  lacking  in  organics  to  interbedded 
silts  and  sands  with  thin  surface  organic  layers.  These  soils  are  generally  well  drained 
(i.e.,  moist)  to  excessively  drained  (i.e.,  dry)  and  exhibit  activities  that  range  from  neutral 
to  slightly  acidic. 

The  mixed  forest  ecotype  is  what  Jorgenson  et  al.,  (2001)  describe  as  Riverine 
Moist  Mixed  Forest.  Vegetation  consists  of  a  range  of  intermediate  successional 
broadleaf  and  needleleaf  species.  Closed  canopies  are  composed  of  Populus 
balsamifera,  Picea  glauca,  and/or  Betula  papyrifera.  Understories  consist  of 
Shepherdia  canadensis,  Dryas  drummondii,  Fragaria  virginiana,  Ceratodon  purpureus, 
Alnus  tenuifolia,  Rosa  acicularis,  Geocaulon  lividum,  Linnaea  borealis  and  Hylocomium 
splendens.  Mixed  ecotype  soils  are  similar  to  those  of  the  broadleaf  ecotype  described 
above. 


The  needleleaf  forest  ecotype  is  what  Jorgenson  et  al.,  (2001)  describe  as 
Riverine  Gravelly  Needleleaf  Forest  and  Moist  Needleleaf  Forest.  Vegetation  consists 
of  a  range  of  late  successional  needleleaf  species.  Open  to  closed  canopies  is 
composed  of  Picea  glauca.  Understories  consist  of  Shepherdia  canadensis,  Solidago 
canadensis,  Geocaulon  lividum  and  Hylocomium  splendens,  Rosa  acicularis.  Ledum 
groenlandicum  and  Calamagrostis  canadensis.  Soil  profiles  exhibit  interbedded  gravel, 
sand  and  silt  layers,  with  thin  surface  organics.  These  soils  are  generally  well  drained 
(i.e.,  moist)  to  excessively  drained  (i.e.,  dry)  and  exhibit  activities  that  are  slightly  acidic. 
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Research  Hypothesis  and  Tests 


The  primary  objective  of  this  study  is  to  describe  and  relate  important  soil 
quality  response/capacity  attributes  that  directly  influence  the  ability  of  the  ecosystem  to 
influence  the  fate  and  transport  of  potential  UXO  contamination.  Given  this  objective, 
this  study  evaluates  ecotype  soil  properties  in  undisturbed  and  disturbed  setting.  Stated 
in  terms  of  hypotheses,  they  are: 

1.  Ho:  There  is  no  relationship  between  ecotype  and  soil  quality. 

Ha:  There  is  a  relationship  between  ecotype  and  soil  quality. 

Test:  Soil  samples  among  ecotypes  and  t-Tests  on  paired  means  of 
selected  indices. 

2.  Ho:  There  is  no  relationship  between  munitions  use  on  a  target  and 
soil  quality  in  proximity  to  the  target. 

Ha:  There  is  a  relationship  between  munitions  use  on  a  target  and  soil 
quality  in  proximity  to  the  target. 

Test:  Soil  samples  around  target  in  concentric  rings  and  t-Tests  on 
paired  means  of  selected  indices. 


Study  Design  and  Execution 

This  design  and  execution  of  this  evaluation  consisted  of  three  phases: 

1 .  Identification  of  study  plots  in  the  Jarvis  Creek  control  area  using  a  land  cover 
classification. 

2.  Soil  sample  collection  in  the  Jarvis  Creek  control  area,  the  Washington  Range 
control  area  and  the  Washington  impact  study  area. 

3.  Soil  sample  processing,  laboratory  and  statistical  analysis. 

Phase  1 :  Land  Cover  Classification 


Initial  efforts  centered  on  creating  a  thematic  land  cover  classification  of  the 
control  and  study  sites  producing  thematic  maps  of  the  Tanana  watershed  in  the  vicinity 
of  Fort  Greely,  Delta  Creek  impact  area,  Washington  Range  impact  area  and  the  Jarvis 
Creek  control  area.  The  objective  was  to  generate  a  land  cover  classification  to  achieve 
the  highest  level  of  accuracy  for  the  general  informational  categories  of  scrub,  broadleaf 
forested,  needleleaf  forested,  mixed  forested  ecotypes  associated  with  these  areas. 
Study  plots  in  the  Jarvis  Creek  control  area  control  site  were  selected  using  these 
maps. 


Landsat  7  ETM+  images  were  reviewed  using  the  USGS  Earth  Explorer  web  site 
(http://edcsns1 7.cr.usas.qov)  using  Latitude  and  Longitude  coordinates  (NW  corner  - 
147.04.00W,  64.20.00N  and  SE  corner  -  145.23.00W,  63.39.00N)  to  obtain  a  listing  of 
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metadata  records  for  all  available  imagery  of  the  area.  Candidate  images  were 
reviewed  based  on  the  lowest  available  aggregate  cloud  cover  and  acquisition  dates 
from  June  to  early  September  to  maximize  temporal  variability  of  vegetation  growth  for 
as  many  of  the  potentially  present  vegetation  classes  while  at  the  same  time  avoiding 
snow  cover  at  higher  elevations.  Additionally,  the  image  had  to  be  recent  acquired 
(i.e.,  no  more  than  2  years  old)  in  order  to  maintain  a  high  degree  of  mapping 
confidence.  A  Landsat  7  ETM+  scene  of  the  area  taken  on  10  September  1999  (path 
67,  row  15)  with  zero  aggregate  cloud  cover  was  obtained  through  US  Army  Space 
Command  (ARSPACE)  Remote  Sensing  unit  in  Colorado  Springs,  Colorado  and  the 
Topographic  Engineering  Center,  Fort  Belvoir,  Maryland. 

Reference  data  for  Fort  Greely,  Alaska  was  obtained  from  the  US  Army  Corps  of 
Engineers  Cold  Regions  Research  and  Engineering  Laboratory  (CRREL),  in  Hanover, 
New  Hampshire.  They  provided  a  draft  report  titled,  “An  Ecological  Land  Survey  for 
Fort  Greely,  Alaska”  authored  by  Jorgenson  et  al.  (2000).  Their  field  data  was  collected 
between  1996  and  1998  and  consists  of  253  surveyed  reference  points  and  126 
surveyed  verification  points  that  describe  the  major  cover  types  found  on  the  Fort 
Greely  reservation. 

The  image  file  was  copied  from  the  CDROMs  to  a  working  directory.  Using 
ERDAS  Imagine  8.4,  the  image  was  display  in  viewer  window.  Header  records,  file 
statistics  and  histograms  showed  that  the  image  was  properly  georectified.  The  scene 
was  subsetted  to  capture  the  area  of  interest  the  Fort  Greely  reservation.  At  this  point, 
no  enhancement  was  made  to  the  image  to  preserve  the  spectral  reflectance 
characteristics  of  the  pixels.  The  scene  had  cloud  cover  and  cloud  shadow  in  the 
southeastern  portion  of  the  image  centered  over  the  Washington  impact  area  (Appendix 
A). 


A  supervised  classification  of  the  scene  was  made  using  the  reference  and 
verification  data  generated  by  CRREL  (Jorgenson  et  al.,  2000).  Reference  point  data 
was  converted  into  an  Arcview  cover,  which  was  then  brought  cover  into  the  Imagine 
viewer.  Attribute  information  was  obtained  for  each  reference  point  using  the  same 
cover  in  an  Arcview  3.2  environment.  Using  the  Imagine  region-growing  tool  for  all  seed 
pixels  (Area  constrained  =  100  pixels  and  Euclidean  Distance  =  20),  pixels  were 
selected  that  corresponded  to  the  appropriate  cover  type.  These  seed  pixels  were  then 
used  to  grow  training  polygons  of  similar  spectral  reflectance.  This  procedure  was 
repeated  for  all  of  the  usable  reference  points.  For  significant  cover  types  without 
reference  points  (e.g.  bare  exposed  rock,  glacier,  snow,  road,  agriculture)  training 
polygons  were  delineated  based  on  a  manual  image  interpretation  of  CIR  composites  of 
the  September  scene. 
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Informational  categories  are  listed  below: 


Level  1  Classification 

Level  2  Classification 

Water 

• 

River 

• 

Pond  and  Lake 

Forest 

• 

Broadleaf 

• 

Needleleaf 

• 

Mixed 

Scrub 

.  Scrub 

Gravel/Pavement 

• 

Roads 

• 

Bare  Exposed  Rock 

• 

Barrens 

Agricultural  Land 

.  Agricultural 

During  this  process,  there  was  some  overlap  of  training  polygons  of  different 
reference  classes.  A  best  estimate  of  which  polygon  to  use  based  upon  a  seed  pixel’s 
proximity  to  an  edge  signature.  A  polygon  whose  seed  pixel  was  best  encompassed  by 
similar  pixels  (i.e.,  farthest  from  pixels  of  different  color)  was  retained.  Reference  points 
associated  with  cloud  and  cloud  shadows  were  not  used.  Upon  completion  of  the 
training  polygon  selection  process  68  polygons  generated  by  seed  pixels  had  been 
obtained  along  with  12  polygons  generated  by  manual  polygon  selection  for  a  total  of  80 
spectral  classes.  A  maximum  likelihood  algorithm  was  used  to  classify  the  whole  image 
based  on  the  80  training  spectral  classes.  A  signature  seperability  analysis  was 
performed  on  the  spectral  classes  using  the  transform  divergence  method  for  all  1-6 
bands  of  information.  This  generated  2,926  pair-wise  comparisons.  Using  the  best 
average  seperability  and  a  cut  off  score  of  1 ,200  (out  of  a  possible  2,000)  spectral 
classes  that  were  similar  were  merged  and  poor  signatures  deleted  bringing  the  training 
signatures  down  to  43  spectral  classes.  Like  cover  classes  (e.g.  Scrub  1  +  Scrub  2  +  ... 
Scrub  /,  etc)  were  then  combined  based  on  histogram  evaluations  using  a  recode 
process  resulting  in  mapping  set  of  21  cover  classes.  Finally,  these  cover  classes  were 
recolored  to  better  distinguish  patterns.  The  classified  image  was  processed  with  a  5X5- 
majority  filter,  which  greatly  enhanced  its  appearance  (Figure  4). 

From  the  maps  of  sampling  areas  of  interest,  over  100  one-hectare  size  discrete 
(i.e.,  non-overlapping  and  homogenous)  sampling  plots,  including  15  per  ecotype,  were 
identified  using  a  random  point  generator  within  ERDAS  Imagine  8.4.  These  sampling 
plots  doubled  as  verification  points.  The  recommended  number  of  verification  points  for 
large  area  classification  should  be  approximately  50  per  ecotype  (Jenson,  1996). 
However,  this  represents  an  enormous  sampling  effort  and  may  be  significantly  more 
than  required  due  to  the  relatively  small  sampling  area  of  the  Jarvis  Creek  control  area. 
Jenson  (1996)  suggests  that  the  total  number  of  verification  points  can  be  estimated  by 
the  formula  for  the  binomial  theory: 

A/=2^  fp)  (a)  where, 

Zis  the  standard  normal  deviate  for  the  95%  two-sided  confidence  interval. 
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p  is  the  expected  percent  error, 
q  =  1 00-p  and 
E  is  the  allowable  error. 

For  this  study,  Z=  1 .96,  p  =85%,  q  =  15%,  and  E  =  10%  yielding  an  A/ =  49 
points.  This  minimum  number  then  would  be  allocated  across  the  ecotypes  of  interest 
as  well  as  any  other  cover  type  of  mapping  interest  (e.g.,  open  or  barren  areas).  Using 
three  ecotypes  (i.e.,  scrub,  broadleaf,  and  mixed)  and  a  cover  type  representing  gravel 
barrens  and  gravel  roads,  this  equates  to  approximately  12  verification  points  per  cover 
type  which  as  mentioned  above  are  the  sampling  plots.  These  sample  plots  were 
circular  with  a  radius  of  55  m  for  an  area  of  one-hectare. 

Phase  2:  Soil  Sample  Collection 

Central  to  the  sampling  effort  is  the  recognition  that  measurements  of 
environmental  phenomena  such  as  soils  or  the  presence  or  absence  of  a  contaminant 
are  generally  characterized  by  great  variability.  In  an  attempt  to  limit  variability  and  the 
associated  uncertainty,  as  well  as  reduce  the  number  of  required  samples,  CRREL  has 
developed  a  composite  sampling  protocol  (Ramsey,  2001,  personal  communication). 
This  method  involves  collecting  numerous  (i.e.,  fifty)  incremental  samples  from  a  given 
area,  compositing  these  incremental  samples  into  one  well  mixed  sample,  then  using 
this  well  mixed  sample  for  all  laboratory  analysis  to  achieve  a  substantial  reduction  in 
the  realized  variability  of  the  data  set.  Since  contamination  is  relatively  more  variable 
than  measures  of  soil  component  characteristics,  a  composting  sampling  approach 
should  work  well  for  soil  property  and  subsequent  soil  quality  assessment.  Therefore, 
all  field  sampling  for  this  research  project  employed  this  approach.  Field  sampling 
consisted  of  two  investigations: 

Investigation  1  -  Control  Area  Characterizations: 

The  purpose  of  this  investigation  is  to  determine  baseline  levels  of  soil  quality 
distribution  related  to  undisturbed  ecotypes  (Scrub,  Broadleaf  Forest,  and  Mixed  Forest) 
referred  to  in  Research  Hypotheses  1 .  Needleleaf  Forest  ecotype  was  not  sampled  due 
to  limited  field  time,  however  nine  plots  was  reconnoitered  for  verification  purposes.  As 
described  above,  the  thematic  classification  for  the  control  area  along  with  a  random 
start  point  generator  was  used  to  determine  discrete  sampling  plots  in  the  control  site. 
Sampling  plots  were  located  and  marked  using  a  Garmin  lll-i-  GPS  receiver.  Each 
sample  consisted  of  fifty  randomly  collected  40  gm  incremental  samples  that  were 
composited  into  one  2  kg  soil  sample  for  that  particular  plot.  Samples  were  collected 
from  the  interface  of  the  organic  and  mineral  stratums  of  the  soil  profile  (i.e.,  the  “O”  and 
“A”  layers  respectively)  within  a  55-meter  radius  of  the  marked  point  (Figures  5a  and  b). 
These  samples  were  air  dried,  sieved  through  a  4  mm  mesh,  mixed  well,  double  bagged 
in  plastic  bags  and  submitted  to  the  CSU  Soils  Laboratory.  At  each  sampling  plot  the 
depth  of  the  cumulative  organic  matter  was  determined  by  digging  a  small  trench  to  the 
depth  of  the  lowest  organic  horizon  or  40  cm  which  ever  was  deeper.  Additionally, 
elevation,  slope  and  aspect,  and  vegetative  cover  descriptions  were  noted  for  each  plot. 
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Investigation  2  -  Physical  Disturbance  in  Washington  Impact  Area: 

The  purpose  of  this  investigation  is  to  determine  soil  quality  in  disturbed  ecotypes 
and  to  characterize  soil  quality  distribution  around  impact  craters  referred  to  in 
Research  Hypotheses  2. 

The  sampling  area  around  an  impact  crater  consisted  of  five  concentric  rings 
radiating  incrementally  outward  from  the  center  of  the  target  (Figures  6a  and  b).  The 
incremental  distance  (/)  from  the  target  array  will  be  determined  by  dividing  the  diameter 
(D)  of  the  observed  disturbed  area  by  five  (i.e.,  /  =  D/5)  creating  five  concentric  rings. 
These  four  points  are  the  sampling  locations  for  each  ring.  Within  each  of  these  rings, 
twenty  40  gm  incremental  samples  were  collected  from  the  top  2  cm  of  the  soil  surface. 
These  4  x  15  =  60  sub-samples  will  be  composited  into  one  2.4  kg  soil  sample  for  that 
particular  sample  ring.  These  samples  were  air  dried,  sieved  through  a  4  mm  mesh, 
mixed  well,  double  bagged  in  plastic  bags  and  submitted  to  the  CSU  Soils  Laboratory. 
Finally,  elevation,  slope  and  aspect  will  be  noted  for  each  target  array. 

Phase  3:  Soil  Sample  Processing.  Laboratory  and  Statistical  Analysis 

Standard  methods  were  used  to  analyzed  soil  samples  and  are  described  in 
Table  3.  Soil  data  and  descriptive  statistics  is  at  Appendix  C. 

Investigation  1  -  Data  and  Analysis 

Twelve  one-hectare  composite  samples  were  taken  from  each  of  the  four 
undisturbed  ecotypes  (i.e.,  Washington  Control  are  scrub  ecotype  and  Jarvis  Creek 
control  area  scrub,  broadleaf,  and  mixed  forest  ecotypes).  A  series  of  random  duplicate 
samples  were  evaluated  to  determine  if  there  were  significant  differences  in  both 
sampling  method  and  laboratory  analysis.  Statistical  analysis  focused  on  detecting 
significant  differences  between  mean  %  organic  matter  and  pH  of  paired  ecotypes. 

Tests  were  performed  using  two-sample  t-Tests  assuming  unequal  sample  population 
variances  and  utilizing  an  alpha  =  0.5. 

Investigation  2  -  Data  and  Analysis 

Data  and  descriptive  statistics  for  the  crater  ring  characterizations  can  be  found 
at  Appendix  D.  Six  craters  were  sampled.  Replicates  were  visually  evaluated  based 
on  error  bars  derived  from  a  95%  confidence  interval  of  their  means.  Statistical  analysis 
of  the  rings  focused  on  detecting  significant  differences  between  mean  %  organic 
matter  of  paired  rings.  Tests  were  performed  using  two-sample  t-Tests  assuming 
unequal  sample  population  variances  and  utilizing  an  alpha  =  0.5. 
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Results  and  Discussion 


Land  Cover  Classification  RS/GIS  Accuracy 

Based  on  vegetative  growth  patterns,  the  classified  image  that  appeared  to 
satisfactorily  classify  cover  types.  The  classified  image  appeared  to  emulate  the 
transitions  between  river,  barren,  scrub,  mixed  forested,  broadleaf  forested,  and 
needleleaf  forested  areas,  and  it  seemed  to  accurately  classified  lakes  and  ponds  and 
agricultural  areas.  Overall  classification  accuracy  was  78.3%  (Table  4).  However, 
classification  of  three  of  the  four  riverine  ecotypes  was  exceptionally  accurate.  Mapping 
of  Scrub,  Broadleaf  and  Mixed  Forest  ecotypes  resulted  in  1 00%,  90%  and  89%  User’s 
Accuracy  respectively  and  the  Kappa  statistics  describing  classified  verses  reference 
data  were  equally  as  good.  Needleleaf  Forest  ecotype  was  not  accurately  classified 
due  to  confusion  with  the  other  three  ecotypes.  This  was  most  probably  due  to 
verification  interpretation  of  the  dominant  canopy  present  during  sampling.  Of  note  was 
the  classification’s  ability  to  map  scrub  ecotype  (User’s  Accuracy  of  1 00%),  which  is  the 
predominant  ecotype  of  the  Washington  impact  area  and  of  other  impact  areas  on  Fort 
Greely. 

This  classification’s  usefulness  may  only  extend  to  riverine  settings  since  the 
preponderance  of  the  verification  points  fell  within  riverine  ecotypes.  Also,  the 
classification  process  grouped  two  or  three  forms  of  ecotypes  (e.g.,  level  3  classification 
of  three  types  of  scrub  vegetation)  which  could  possibly  have  been  labeled  in  further 
detail  had  additional  information  been  available.  Even  so,  this  method  may  still  be 
appropriate  for  impact  area  soil  quality  evaluations  in  other  settings  as  long  as  larger 
scale  ecosystem  components  such  as  physiography  are  held  constant. 

The  use  of  remotely  sensed  data  in  the  form  of  Landsat  7  ETM+  coupled  with 
computer  aided  analysis  for  classification  into  a  geographic  information  system  appears 
to  be  a  viable  method  of  determining  ecotype  in  a  riverine  setting.  When  soil  quality  is 
correlated  with  ecotype,  it  should  be  a  practical  method  of  inferring  soil  quality  through 
the  use  of  remotely  sensed  data  for  situations  in  which  large  inaccessible  areas  must  be 
evaluated.  Therefore,  it  should  be  of  significant  value  in  determining  soil  quality  for 
similar  investigations. 


Indicators  of  Soil  Quality 

A  summary  of  statistical  tests  and  results  is  at  Table  5.  Since  adsorption  and 
chelation  are  directly  related  to  cation  exchange  capacity  (CEC),  and  CEC  is  a  function 
of  the  colloidal  content  contributed  by  both  organic  matter  and  clays,  a  regression  of 
CEC  on  both  %  organic  matter  and  %  clay  was  performed  to  determine  which  was  the 
dominant  contributor.  CEC  proved  to  be  highly  correlated  to  %  organic  matter  with  77% 
of  the  variability  in  CEC  represented  by  %  organic  matter  within  the  control  areas  and 
with  86%  of  the  variability  in  CEC  represented  by  %  organic  matter  within  the  impact 
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area.  Clay  was  not  correlated  with  CEC.  Therefore,  %  organic  matter  was  used  to 
evaluate  soil  quality.  Additionally,  soil  pH  affects  CEC,  so  it  too  was  used  as  an 
indicator  of  soil  quality. 


Soil  Quality  Characteristics  of  Undisturbed  Ecotypes 

Statistical  analysis  of  the  Jarvis  Creek  control  area  ecotypes  showed  no 
evidence  that  %  organic  matter  or  pH  was  different  among  these  ecotypes  (Table  6). 
However,  there  may  be  soil  quality  variables  that  co-vary  within  ecotypes  that  may  be 
used  to  detect  significant  differences  between  them.  The  Washington  Range  control 
area  scrub  ecotype  %  organic  matter  and  pH  were  significantly  different  from  each  of 
the  Jarvis  Creek  control  area  ecotypes.  These  results  indicate  that  there  may  be  a 
significant  difference  in  the  soil  quality  of  the  Jarvis  Creek  control  area  ecotypes  and 
Washington  Range  control  area  scrub  ecotype.  This  may  have  implications  for  the  use 
of  remotely  sensed  ecotype  spectral  signatures  in  assessing  soil  quality,  but  the  results 
may  have  been  better  if  the  Washington  range  and  impact  areas  were  included  in  the 
classification  process.  Also,  a  level  3  classification  may  have  resulted  in  better 
discrimination  between  different  vegetative  forms  within  ecotypes  (e.g,  differentiating 
between  a  scrub  ecotype  dominated  by  Dryas  drummondii  and  one  dominated  by 
Populus  baisamifera  saplings).  Soil  quality  differences  between  these  forms  may  then 
be  more  pronounced. 


Soil  Quality  Characteristics  of  Undisturbed  and  Disturbed  Ecotypes  (Craters) 

Comparison  of  undisturbed  Washington  Range  control  area  scrub  to  the 
Washington  Impact  area  scrub  showed  no  significant  differences  in  %  organic  matter 
and  pH.  There  was  a  significant  difference  between  soil  quality  of  craters  and  sampled 
plots  within  the  impact  area  as  well  as  a  significant  difference  between  soil  quality  of 
craters  and  Washington  Range  control  area  scrub.  Statistical  analysis  of  the  crater  ring 
replicates  showed  that  both  sampling  technique  and  laboratory  analysis  appear  to  be 
consistent  and  reasonably  precise  (Figures  7a-f).  Soil  nutrient  values  are  the  most 
variable,  but  key  measures  of  adsorption  capacity  (i.e.,  CEC  and  pH)  are  very 
consistent.  Data  analysis  shows  that  as  one  moves  from  the  center  of  the  crater  to  the 
outer  most  sampled  ring  there  is  a  significant  gradient  of  soil  quality  characteristics 
(Table  7).  While  adjacent  rings  do  not  show  significant  differences  in  %  organic  matter, 
there  is  a  significant  difference  between  rings  1  &  3,  2  &  4,  and  3  &  5.  This  difference  is 
apparent  in  graphs  of  %  organic  matter  verses  distance  from  the  crater  center  (Figure  8 
a-f). 


This  analysis  shows  that  cratering  is  a  significant  physical  disturbance  within  the 
impact  area.  If  the  frequency  and  intensity  of  this  physical  disturbance  exceeds  the 
capacity  of  the  system  to  recover,  the  system  may  reach  a  threshold  from  which  both 
recovery  and  capacity  for  contaminant  attenuation  is  lost.  Whether  the  disturbance 
regime  has  or  will  reach  a  critical  resiliency  threshold  can  be  evaluated  by  comparing 
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crater  age.  Relatively  higher  concentrations  of  organic  matter  were  found  in  the  center 
of  craters  that  were  four  years  or  older.  This  is  a  reasonable  finding  given  the  nature  of 
the  soil  forming  processes  at  work  in  the  Washington  Impact  Area.  Craters  provide  a 
depression  that  can  capture  wind  blown  leaves,  silts  and  finer  organic  particulates,  as 
well  as  seeds.  They  may  also  capture  snow  and  serve  as  a  moisture  source  for  plant 
establishment.  Most  of  the  craters  sampled  had  small  balsam  poplar  saplings  at  their 
center  (used  to  approximate  the  age  of  the  craters).  It  appears  that  the  physical 
disturbance  and  displacement  of  soil  material  eventually  serves  as  a  catalyst  for  organic 
matter  additions  assisting  recovery.  Also,  the  %  organic  matter  of  the  surface  soil 
appears  to  recover  as  shown  by  the  flattening  of  the  gradient  curve  on  older  craters. 


Conclusion 


Follow-on  Analysis 

A  principle  component  analysis  of  soil  properties  may  provide  greater  insight  into 
other  important  ecotype  associated  soil  properties.  While  the  statistical  methods  used 
in  this  analysis  did  not  detect  a  significant  difference  in  the  soil  quality  of  undisturbed 
Jarvis  Creek  control  area  ecotypes  for  %  organic  matter,  a  principle  component  analysis 
may  provide  better  resolution  of  differences.  Also,  the  data  may  be  better  suited  for 
non-parametric  analysis. 

CRREL  gathered  munitions  contamination  data  for  the  Washington  Impact  Area 
during  the  same  period  this  study  gathered  soil  quality  samples.  Many  of  the  CRREL 
samples  were  superimposed  over  soil  quality  sampling.  The  CRREL  data  should  be 
available  in  early  2002.  Again,  utilizing  a  principle  component  analysis  and  non- 
parametric  methods  of  analysis  may  yield  additional  findings.  If  contamination  can  be 
correlated  with  soil  quality,  in  particular  to  an  indicator  such  as  %  organic  matter  which 
appears  to  be  a  reasonably  consistent  metric  across  all  sampling  done  in  this  study, 
then  it  may  be  possible  to  rate  a  given  soil  for  its  resilience  and  attenuation  ability. 


Development  and  Use  of  Soil  Organic  Matter  Resiliency  Metric 

The  soil  component  performs  the  following  functions:  provides  habitat  for  soil 
microorganisms,  serves  as  an  engineering  medium  and  structural  support  for  plants, 
acts  to  filters  water  resources,  and  serves  as  a  nutrient  bank  for  vegetation  (Bradey  and 
Weil,  1999).  Quality  soils  have  been  described  by  their  soil  organic  matter.  Soil  organic 
matter  (SOM)  is  a  principal  element  of  the  soil  component  affecting  all  of  the  above. 
SOM  enhances  microbial  functions  increasing  energy  flux  and  rate  of  nutrient  cycling, 
which  in  turn  aids  degradation  of  organic  compounds  resulting  in  less  soil  and  water 
pollution  and  greater  plant  production.  SOM  is  the  source  of  humic  substances  that  can 
increased  a  soil’s  buffer  capacity  resulting  in  a  more  stable  soil  pH.  SOM  can  increase 
water-holding  capacity,  which  provides  for  greater  plant  production.  SOM  is  the  source 
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of  non-humic  substances  that  can  increase  aggregate  stability,  which  means  less  soil 
erosion  and  less  soil  and  water  pollution.  These  same  non-humic  substances  can 
provide  increased  adsorption  and  metal  chelation  capacity  enabling  the  soil  component 
to  retain  organic  and  inorganic  pollutants,  which  also  contributes  to  less  soil  and  water 
degradation. 

It  is  these  properties,  and  in  particular  SOM’s  ability  to  influence  the  fate  and 
transport  of  contaminants  through  an  ecosystem,  that  provides  the  basis  for  future  work. 
If  one  can  characterize  the  distribution  of  SOM,  predict  its  dynamics,  and  ascribe 
measures  of  its  attenuation  potential,  then  the  implications  for  adaptive  ecosystem 
management  in  the  face  of  evolving  contaminant  threats  may  be  profound.  Resource 
managers  and  environmental  scientists  could  adapt  their  activities  to  maximize  this 
potential  thereby  mitigating  or  avoiding  disturbances  that  are  beyond  an  ecosystem’s 
resistance  and  resilience  capability. 

Agricultural  practices  such  as  tillage,  fertilization,  pesticide  application  and 
cropping  have  been  universally  recognized  as  anthropogenic  disturbances  of  varying 
degrees  affecting  ecosystem  productivity.  SOM  is  central  to  understanding  the  effects  of 
these  practices.  This  agricultural  scenario  can  readily  be  adapted  toward  range  and 
impact  area  ecosystem  management.  Recently,  the  National  Research  Council,  (2000) 
has  proposed  that  SOM  should  be  widely  recognized  as  an  indicator  of  soil  quality  and 
productivity  because  of  its  ability  to: 

.  provide  a  nutrient  and  energy  source  for  soil  biota 
.  improve  soil  structure  by  strengthening  soil  aggregates 
.  increase  water  holding  capacity 

.  promote  infiltration  and  reduce  erosion 

.  increase  cation  exchange  capacity 
.  chelate  metals 

.  Influence  the  fate  and  transport  of  pesticides 

Also,  SOM  is  considered  the  single  most  important  indicator  of  soil  health  (Larson  and 
Pierce,  1994).  SOM  measurable  attributes  are  listed  in  Table  8  (Fenton  et  al,  1999). 
They  provide  a  sensitive  measure  of  SOM’s  contribution  to  an  ecosystem’s  soil  function. 

Means  of  Simulating  Soil  Organic  Matter  Dynamics 

Once  SOM  for  a  particular  system  Is  characterized  it  can  then  be  readily  modeled 
and  used  as  a  resiliency  metric.  The  next  step  in  exploring  ecosystem  resilience  in  the 
face  of  munitions  disturbance  is  to  use  SOM  dynamics  output,  or  the  fluxes  in  SOM 
pools  as  a  result  of  natural  and  anthropogenic  driving  variables,  to  provide  an 
established  metric  for  determining  an  ecosystem’s  response  to  various  management 
practices.  The  findings  of  this  study  should  enable  effective  parameterization  of  such  a 
simulation. 
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Table  1.  Notional  Range  Requirements  for  US  Army  Installations  (TC  25-8).  Target  and  target  array  calculations  are  based  on  ranges  descriptions  in  TC  25-8. 
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Table  2.  Summary  of  area  determinations  and  potential  contaminant  loading  rates. 
Tabular  data  can  be  found  at  Appendix  B. 


Acre 

Hectare 

Ft  Greely  Installation  Area 

662000 

267905 

Total  High  Hazard  Impact  Area 

85042 

34416 

Surface  Area  of  Target  Arrays 

2409.3 

975 

Average  Area  Cratered  per  Year 

25.324 

10.25 

Percent  Area 

%  Ft  Greely  that  is  High  Hazard  Impact  Area 

12.85% 

%  High  Hazard  Impact  Area  in  Use 

2.83  % 

%  Target  Array  Surface  Area  Cratered  per  Year 

1.05% 

%  High  Hazard  Impact  Area  Cratered  per  Year 

0.0298  % 

Loading  Rate 

Total  Mass  of  Explosive  Filler  (DUD  +  LO) 

4460  kg 

Theo.  Max.  Explosive  Loading  Rate  TGT  Arrays 

0.46  g/m2/yr 

Theo.  Max.  Explosive  Loading  Rate  Crater  Area 

43.52  g/m2/yr 

Soil  properties: 

Methods: 

Reference: 

pH 

Soil  to  water  ratio  of  1 :1 

2,  p.487 

pH  Buffer 

SMP  Buffer 

2,  pp.502-505 

Conductivity 

Soil  to  water  ratio  of  1 :1 

2,  pp.420-427 

Cation  Exchange 

Measure  Na;  account  for 
entrained  Na  with  nitrite. 

2,  pp.  1201 -1230 

Texture 

Hydrometer 

1 ,  method  1 5-4 

%  Organic  Carbon 

Modified  Walkely-Black 

2,  pp.995-996 

%  Organic  Matter 

1 .724  X  %  Organic  Carbon 

Soil  Nutrients: 

Nitrogen  as  ammonia 
Nitrogen  as  nitrate 

2M  potassium  chloride 
extract,  automated 
phenate  method  by  FIA 

2,  pp.1 146-1 162 

Total  N 

CHN  furnace 

2,  pp.1 085-1 122 

Total  C 

CHN  furnace 

2,  pp.967-977 

Phosphorous 

Mehlich  3  Extract 

2,  pp.893-894 

Potassium 

Mehlich  3  Extract 

2,  pp.893-894 

Calcium 

Mehlich  3  Extract 

2,  pp.893-894 

Magnesium 

Mehlich  3  Extract 

2,  pp.893-894 

Extractable  Metals: 

Cadmium 

Mehlich  3  Extract 

2,  pp.893-894 

Chromium 

Mehlich  3  Extract 

2,  pp.893-894 

Copper 

Mehlich  3  Extract 

2,  pp.893-894 

Lead 

Mehlich  3  Extract 

2,  pp.893-894 

Zinc 

Mehlich  3  Extract 

2,  pp.893-894 

Table  3.  List  of  soil  properties,  nutrients  and  metals  analysis,  and  laboratory 
methods. 


Reference: 

1 .  Klute,  A.  1986.  Methods  of  Soil  Analysis:  Part  1 .  Physical  and  Mineralogical 
Methods.  Soil  Science  Society  of  America,  Inc.,  Madison,  Wisconsin. 

2.  Spark,  D.L.,  (ed).  1996.  Methods  of  Soil  Analysis:  Part  3.  Chemical  Methods. 
Soil  Science  Society  of  America,  Inc.,  Madison,  Wisconsin. 


CLASSIFICATION  ACCURACY  ASSESSMENT  REPORT 


ERROR  MATRIX 


Reference  Data  1 

Classified  Data 

Mixed 

Forest 

Scrub 

Needleleaf 

Forest 

Gravel 

Pavement 

Broadleaf 

Forest 

Mixed  Forest 

8 

0 

1 

0 

0 

Scrub 

0 

14 

0 

0 

0 

Needleleaf  Forest 

3 

4 

7 

0 

2 

Gravel/Pavement 

0 

0 

1 

9 

1 

Broadleaf  Forest 

1 

0 

0 

0 

9 

Column  Total 

12 

18 

9 

9 

12 

ACCURACY  TOTALS 


Class  Name 

Reference 

Classified 

Number 

Producers 

Users 

Totals 

Totals 

Correct 

Accuracy 

Accuracy 

Mixed  Forest 

12 

9 

8 

66.67% 

88.89% 

Scrub 

18 

14 

14 

77.78% 

100.00% 

Needleleaf  Forest 

9 

16 

7 

77.78% 

43.75% 

Gravel/Pavement 

9 

11 

9 

100.00% 

81.82% 

Broadleaf  Fores 

12 

10 

9 

75.00% 

90.00% 

Totals  or  Average 

60 

60 

47 

79.45% 

80.89% 

Overall  Classification  Accuracy  =  78.33% 


KAPPA  (K'^)  STATiSTiCS 
Overaii  Kappa  Statistics  =  0.7289 

Conditional  Kappa  for  each  Category. 


Class  Name 

Kappa 

Mixed  Forest 

0.8611 

Scrub 

1 

Needleleaf  Forest 

0.3382 

Gravel/Pavement 

0.7861 

Broadleaf  Forest 

0.875 

Table  4.  Error  matrix,  Mapping  Accuracy  and  Kappa  Statistics  for  Ecotype  Classification 
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Table  6.  Undisturbed  ecotype  analysis.  Shaded  blocks  depict  ecotypes  for  which  there 
is  no  evidence  of  a  difference  between  ecotype  for  %  organic  matter.  “Scrub  W”  is 
Washington  Range  Control  Area  scrub  and  “Scrub  J”  is  Jarvis  Creek  Control  Area 
Scrub. 


Table  7.  Crater  ring  analysis  in  which  common  crater  rings  were  averaged  together. 
Shaded  blocks  depict  rings  for  which  there  is  no  evidence  of  a  difference  between  rings 
for  %  organic  matter. 


Table  8.  SOM  attributes,  measures  and  environmental  effects  summarized  from  (Fenton 
et  al.,  1999). 


SOM  Attribute 

Measure(s) 

Environmental  Effect 

Buffer  Capacity 

Buffer  Capacity 

Stable  pH  to  aid  adsorption 

Water  Holding 

Capacity 

Bulk  Density  and  Texture 

Biomass  production 

Aggregate  Stability 

Bulk  Density  and  Texture 

Erosion  Prevention 

Adsorption  Capacity 

Cation  Exchange  Capacity  and  pH 

Organic  and  Inorganic 

Adsorption 

Study  and  Control  Areas 


Figure  1 .  Map  of  the  Fort  Greely  Reservation  depicting  the  Jarvis  Creek  and 
Washington  Range  Control  Areas  as  well  as  the  Washington  Impact  Area. 
Reference  and  verification  points  are  sampling  points  described  by  Jorgenson  et 
al.,  2001. 


Figure  2.  Information  management  integration  of  spatial  scaling.  Management 
resolution  and  Landsat  multispectral  data  comparisons  yield  scale  of  observation  for 
sampling  plots. 


Fort  Greely,  Alaska  Climograph 


“O—Normal  Dailey  Mean  Temperature,  Deg  C 
^>"Normal  Monthly  Precipitation,  mm 


Average  Temperatures  above  10  C,  "Approx.  Growing  Season 
Average  Temperatures  Below  0  C,  “Approx.  Period  Frozen  Soils" 


Figure  3.  Climate  graph  for  Fort  Greely  Alaska  depicting  effective  growing  season 
(Diersing  et  a.,  1990).  Climate  data  obtained  from  Big  Delta  FAA/AMOS  AP,  AK 
(500770).  2001. 
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Figure  4.  Jarvis  Creek  Control  Area  Map  with  Sampling  Plots 


570000 


Sample  Plot  Random  Walk 


•  Mark  plot  center  point 

•  Stay  within  a  radius  of 
55  meters 

•  Randomly  walk  around 
center  point 

•  Collect  fifty  40  g 
incremental  samples  = 
2  kg  sample 

•  Organic  matter  depth 
from  center  point. 


Figure  5a.  Depiction  of  hectare  plot  sampling  method  using  random  walk  collecting  50 
incremental  samples  for  one  composite  sample. 


Figure  5b.  Example  of  ‘O’  and  ‘A’  horizon  interface  from  which  soil  samples  were 
collected. 


Figure  6a.  Depiction  of  crater  sampling  method  using  random  walk  collecting  20 
incremental  samples  for  one  composite  ring  sample. 


Figure  6b.  Actual  layout  of  crater  sampling  during  fieldwork.  This  is  Crater  7  which  is 
estimated  to  be  a  1  to  2  year  old  1 05mm  crater. 


Sample  ID  # 

pH 

Buffer  pH 

CEC  (meq 
/lOOg) 

Total  N 
(ppm) 

P  (ppm) 

K  (ppm) 

Total  % 
N 

Total  % 
C 

CR  9-3-1 

7.6 

7.5 

5.3 

23.8 

9.13 

27.8 

0.581 

CR  9-3-2 

7.8 

7.5 

5.4 

6.2 

10.8 

31.9 

0.827 

CR  9-3-3 

7.5 

7.5 

7.3 

37.3 

9.39 

36.8 

0.044 

0.568 

Confidence  Level (95.0%) 

0.3794586 

0 

2.799483 

38.740099 

2.2742298 

11.283637 

0.0222651 

0.3625 

- 1 

Crater  9,  Ring  3 


Crater  9,  Ring  3 


Analyte 


0  Total  N  (ppm) 
EIP  (ppm) 
m  K  (ppm) 


1 


2 


3 


Error  bars  are  based  on  a  95%  confidence  level  for  the  mean 
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Figure  7a.  Crater  9  Ring  3  analysis  of  replicates. 


Sample  ID  # 

pH 

Buffer  pH 

CEC  (meq 
/100g) 

Total  N 
(ppm) 

P  (ppm) 

K  (ppm) 

Total  % 
N 

Total  % 
C 

CR  2-5-1 

7.4 

7.5 

19.5 

29.1 

9.80 

45.6 

0.157 

1.948 

CR  2-5-2 

7.3 

7.5 

26.2 

31.9 

10,1 

66.9 

0.138 

1.768 

CR  2-5-3 

CO 

4.2 

11.8 

11.1 

52.0 

0.160 

2.458 

Confidence  Level(95.0%) 

0 

28.012826 

27.044434 

1.6620236 

27.214542 

0.0296367 

0.8890998 
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Figure  7c.  Crater  2  Ring  5  analysis  of  replicates. 


Sample  ID  # 

pH 

Buffer  pH 

CEC  (meq 
/lOOg) 

Total  N 
(ppm) 

P  (ppm) 

K  (ppm) 

Total  % 
N 

Total  % 
C 

CR  17-1-1 

7.6 

7.5 

10,1 

21.8 

1.39 

42.2 

0.045 

0.908 

CR  17-1-2 

7.7 

7.5 

11.1 

9.2 

5,12 

92.6 

0.049 

0.828 

CR  17-1-3 

7.5 

7.5 

10.4 

16.8 

4.91 

40.4 

0.050 

0.831 

Confidence  Level (95.0%) 

0.2484139 

0 

1.2747614 

15.760751 

5.2050531 

73.601083 

0.0065724 

0.1126478 
—  r 

Crater  17,  Ring  1 
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Figure  7d.  Crater  17  Ring  1  analysis  of  replicates. 


Sample  ID  # 


CR  16-2-1 


CR  16-2-2 


CR  16-2-3 


CEC  (meq 
/lOOg) 

Total  N 
(ppm) 

P  (ppm) 

K  (ppm) 

Total  % 
N 

Total  % 
C 

4.1 

50.5 

2.62 

71.0 

0.036 

16.3 

32.5 

2.81 

69.3 

0.112 

14.1 

53.7 

3.22 

84.0 

0.078 

1.121 

16.152638 

28.390283 

0.7550735 

20.098377 

0.0945715 

1.4913618 

Crater  7  Organic  Matter 
1 05mm,  1  to  2  years  old 

Figure  8a.  Crater  7  change  in  %  organic  matter  as  a  function  of  distance  from  crater  center. 


Mk82,  4-5  years  old 

Figure  8b.  Crater  1 5  change  in  %  organic  matter  as  a  function  of  distance  from  crater  center 


Landsat  7  Enhanced  Thematic  Mapper  Image 
Color  IR,  Fort  Greely  and  Vicinity 
10  Spetember  1999 
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Appendix  A 


Appendix  B.  Tabular  Data  for  Munitions  Loading 


Washington  Range  Control  Area  Scrub  Ecotype  14-12-01 
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Jarvis  Creek  Control  Area  Mixed  Forest  Ecotype  14-12-01 
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Washington  Impact  Area  Craters  Washington  Scrub  Ecotype  14-12-01 
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Appendix  C 


CEMML  Soil  Data  Regression  14-12-01 

Jarvis  Creek  and  Cation  Exchange  Capacity  on  Organic  Matter 

Washington  Range 
Control  Areas 


SUMMARY  OUTPUT 


1  Repression  Statistics 

Multiple  R 

0.878704464 

R  Square 

0.772121534 

Adjusted  R  Squa 

0.770670079 

Standard  Error 

6.416563842 

Observations 

159 

ANOVA 

df 

SS 

MS 

F 

Sipnificance  F 

Regression 

1 

21902.1662 

21902.1662 

531.9637403 

2.74386E-52 

Residual 

157 

6464.049772 

41.17229154 

Total 

158 

28366.21597 

PARAMETER  ESTIMATES 


Coefficients 

standard  Error 

tStat 

P-value 

Lower  95% 

4.98805884 

0.688799478 

7.241670465 

1.86558E-11 

3.627548604 

6.348569075 

|x  Variable  OM 

2.422021342 

0.105011521 

23.06433915 

2.74386E-52 

2.214603576 

2.629439109 

RESIDUAL  OUTI 

=UT 

Observation 

Plot 

Observed  Y 

Predicted  Y 

Residuals 

Standard 

Residuals 

1 

108 

2.2 

6.704 

7.596 

1.188 

2 

2501 

2.4 

16.103 

2.440 

0.382 

3 

CR7-1-3 

3.1 

6.546 

-3.538 

-0.553 

4 

CR7-1-2 

3.2 

6.763 

-2.468 

-0.386 

5 

CR  17-3-1 

3.4 

6.445 

-2.268 

-0.355 

6 

CR  17-3-3 

3.5 

6.345 

-3.439 

-0.538 

7 

CR7-2-3 

3.9 

6.441 

-3.611 

-0.565 

8 

CR7-1  -1 

4.0 

7.410 

-2.090 

-0.327 

9 

CR1 5-3-3 

4.1 

6.813 

-3.270 

-0.51 1 

10 

CR  15-4-2 

4.1 

6.821 

-4.504 

-0.704 

11 

CR  16-2-1 

4.1 

6.658 

-6.675 

-1.044 

12 

CR7-3-1 

4.1 

6.871 

-4.270 

-0.668 

13 

CR2-5-3 

4.2 

15.252 

-3.300 

-0.516 

14 

CR7-2-1 

4.2 

6.779 

-3.904 

-0.610 

15 

CR9-2-2 

4.2 

7.393 

-3.861 

-0.604 

16 

CR1  5-2-2 

4.3 

7.548 

-4.647 

-0.727 

17 

CR  15-3-1 

4.4 

6.905 

-3.694 

-0.577 

18 

CR7-2-2 

4.5 

7.059 

-0.807 

-0.126 

19 

E07 

4,9 

5.243 

3.552 

0.555 

20 

CR15-1-3 

4.9 

9.886 

5.757 

0.900 

21 

CR  15-4-3 

4.9 

7.435 

2.047 

0.320 

22 

CR  17-4-1 

4.9 

7.230 

9.203 

1.439 

23 

CR  15-4-1 

5.1 

7.452 

4.692 

0.733 

24 

CR  15-3-2 

5.2 

8.103 

9.515 

1.488 

25 

CR9-3-1 

5.3 

7.414 

-0.369 

-0.058 

26 

CR  15-2-3 

5,4 

8.036 

1.944 

0.304 

27 

CR7-3-3 

5,4 

7.468 

2.817 

0.440 

28 

CR9-3-2 

5.4 

8.441 

2.538 

0.397 

Filename:  soil  data.xls 
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t-Tests 

Comparing  %  Organic  Matter  Means 


14-12-01 


CEMML  Soil  Data 
Jarvis  Creek  and 
Washington  Range 
Control  Areas 


Ecotype 

Scrub  W 

Scrub  J 

Broadleaf 

Mixed 

Code 

A 

B 

C 

D 

Comparison 

A-B 

- 

- 

A-C 

B-C 

-- 

- 

A-D 

B-D 

C-D 

- 

Ecotype 

Scrub  W 

Scrub  J 

Broadleaf 

Mixed 

4.98 

2.01 

8.68 

5.58 

3.28 

2.63 

7.66 

19.10 

3.34 

2.73 

9.93 

11.36 

3.38 

1.95 

17.26 

24.95 

1.76 

4.06 

14.81 

6.63 

%  Organic  Matter 

3.73 

18.86 

15.19 

5.48 

2.73 

8.83 

11.18 

12.01 

2.10 

17.93 

10.21 

16.52 

0.71 

13.13 

7.44 

10.64 

4.78 

13.89 

7.51 

7.89 

2.10 

15.46 

14.67 

11.54 

3.10 

23.26 

10.05 

5.97 

Descriptive  Statistics 

Scrub  W 

Scrub  J 

Broadleaf 

Mixed 

Mean 

2.998036 

10.394 

11.21678 

11.47164 

Standard  Error 

0.35349 

2.206341 

0.985813 

1.74832 

Median 

3.191986 

10.97929 

10.13022 

10.99912 

Mode 

2.098108 

#N/A 

#N/A 

#N/A 

Standard  Deviation 

1 .224524 

7.642989 

3.414958 

6.056357 

Sample  Variance 

1.499458 

58.41528 

11.66194 

36.67946 

Kurtosis 

-0.02208 

-1.47804 

-1.14837 

0.784335 

Skewness 

-0.08112 

0.241315 

0.549006 

1.116313 

Range 

4.27552 

21.30692 

9.816456 

19.46913 

Minimum 

0.708564 

1.949844 

7.440784 

5.477148 

Maximum 

4.984084 

23.25676 

17.25724 

24.94628 

Sum 

35.97643 

124.728 

134.6013 

137.6597 

Count 

12 

12 

12 

12 

Conf  Level(95.0%) 

0.778026 

4.856126 

2.169762 

3.848027 

Summary  of  t-Tests 

Scrub  W 

Scrub  J  Broadleaf  Mixed 

No  evidence  to  reject  Ho: 

Means  are  equal  I  I 


Filename:  undisturbed  ecotypes.xls 
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CEMML  Soil  Data 
Jarvis  Creek  and 
Washington  Range 
Control  Areas 


t-Tests 

Comparing  %  Organic  Matter  Means 


14-12-01 


B  to  C  %  OM 


Scrub  J 

Broadleaf 

2.01 

8.68 

2.63 

7.66 

2.73 

9.93 

1.95 

17.26 

4.06 

14.81 

18.86 

15.19 

8.83 

11.18 

17.93 

10.21 

13.13 

7.44 

13.89 

7.51 

15.46 

14.67 

23.26 

10.05 

t-Test:  Two-Sample  Assuming  Unequal  Variances 


Scrub  J  Broadleaf  \ 

Mean 

10.394 

11.21678 

Variance 

58.41528 

11.66194 

Observations 

12 

12 

Hypothesized  Mean  Difference 

0 

df 

15 

tStat 

-0.34048 

P(T<=t)  one-tail 

0.3691 1 

t  Critical  one-tail 

1.753051 

P(T<=t)  two-tail 

0.73822 

t  Critical  two-tail 

2.131451 

Ho:  Means  are  equal  (difference  =  0) 

Ha:  Means  are  not  equal  (differnce  not  equal  to  0) 
alpha  =  .05 

Conclusion:  No  evidence  to  reject  Ho 


Filename:  undisturbed  ecotypes.xls 
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CEMML  Soil  Data 
Jarvis  Creek  and 
Washington  Range 
Control  Areas 


t-Tests 

Comparing  %  Organic  Matter  Means 


14-12-01 


A  to  B  %  OM 


Scrub  W 

Scrub  J 

4.98 

2.01  ■ 

3.28 

2.63 

3.34 

2.73 

3.38 

1.95 

1.76 

4.06 

3.73 

18.86 

2.73 

8.83 

2.10 

17.93 

0.71 

13.13 

4.78 

13.89 

2.10 

15.46 

3.10 

23.26 

t-Test:  Two-Sample  Assuming  Unequal  Variances 


Scrub  W 

Scrub  J 

Mean 

2.998036 

10.394 

Variance 

1.499458 

58.41528 

Observations 

12 

12 

Hypothesized  Mean  Difference 

0 

df 

12 

tStat 

-3.30993 

P(T<=t)  one-tail 

0.003113 

t  Critical  one-tail 

1.782287 

P(T<=t)  two-tail 

0.006225 

t  Critical  two-tail 

2.178813 

Ho:  Means  are  equal  (difference  =  0) 

Ha:  Means  are  not  equal  (differnce  not  equal  to  0) 
alpha  =  .05 

Conclusion:  Reject  Ho,  means  are  not  equal. 


Filename:  undisturbed  ecotypes.xls 
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CEMML  Soil  Data 
Jarvis  Creek  and 
Washington  Range 
Control  Areas 


t-Tests 

Comparing  %  Organic  Matter  Means 


14-12-01 


B  to  D  %  OM 


Scrub  J 

Mixed 

2.01 

5.58 

2.63 

19.10 

2.73 

11.36 

1.95 

24.95 

4.06 

6.63 

18.86 

5.48 

8.83 

12.01 

17.93 

16.52 

13.13 

10.64 

13.89 

7.89 

15.46 

11.54 

23.26 

5.97 

t-Test:  Two-Sample  Assuming  Unequal  Variances 


Scrub  J  Mixed  \ 

Mean 

10.394 

11.47164 

Variance 

58.41528 

36.67946 

Observations 

12 

12 

Hypothesized  Mean  Difference 

0 

df 

21 

tStat 

-0.38281 

P(T<=t)  one-tail 

0.352853 

t  Critical  one-tail 

1.720744 

P(T<=t)  two-tail 

0.705705 

t  Critical  two-tail 

2.079614 

Ho:  Means  are  equal  (difference  =  0) 

Ha:  Means  are  not  equal  (differnce  not  equal  to  0) 
alpha  =  .05 

Conclusion:  No  evidence  to  reject  Ho 


Filename:  undisturbed  ecotypes.xls 
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CEMML  Soil  Data 
Jarvis  Creek  and 
Washington  Range 
Control  Areas 


t-Tests 

Comparing  %  Organic  Matter  Means 


A  to  C  %  OM 


Scrub  W 

Broadleaf 

4.98 

8.68 

3.28 

7.66 

3.34 

9.93 

3.38 

17.26 

1.76 

14.81 

3.73 

15.19 

2.73 

11.18 

2.10 

10.21 

0.71 

7.44 

4.78 

7.51 

2.10 

14.67 

3.10 

10.05 

t-Test:  Two-Sample  Assuming  Unequal  Variances 


Scrub  W 

Broadleaf 

Mean 

2.998036 

11.21678 

Variance 

1.499458 

11.66194 

Observations 

12 

12 

Hypothesized  Mean  Difference 

0 

df 

14 

t  Stat 

-7.84774 

P(T<=t)  one-tail 

8.56E-07 

t  Critical  one-tail 

1.761309 

P(T<=t)  two-tail 

1.71E-06 

t  Critical  two-tail 

2.144789 

Ho:  Means  are  equal  (difference  =  0) 

Ha:  Means  are  not  equal  (differnce  not  equal  to  0) 
alpha  =  .05 

Conclusion:  Reject  Ho,  means  are  not  equal. 


14-12-01 


Filename:  undisturbed  ecotypes.xls 
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CEMML  Soil  Data 
Jarvis  Creek  and 
Washington  Range 
Control  Areas 


t-Tests 

Comparing  %  Organic  Matter  Means 


C  to  D  %  OM 


Broadleaf 

Mixed 

8.68 

5.58 

7.66 

19.10 

9.93 

11.36 

17.26 

24.95 

6.63 

15.19 

5.48 

11.18 

12.01 

10.21 

16.52 

10.64 

7.51 

7.89 

14.67 

11.54 

10.05 

5.97 

t-Test:  Two-Sample  Assuming  Unequal  Variances 


Broadleaf  Mixed 

Mean 

11.21678  11.47164 

Variance 

11.66194  36.67946 

Observations 

12  12 

Hypothesized  Mean  Difference 

0 

df 

17 

tStat 

-0.12698 

P(T<=t)  one-tail 

0.450222 

t  Critical  one-tail 

1.739606 

P(T<=t)  two-tail 

0.900445 

t  Critical  two-tail 

2.109819 

Ho:  Means  are  equal  (difference  =  0) 

Ha:  Means  are  not  equal  (differnce  not  equal  to  0) 
alpha  =  .05 

Conclusion:  No  evidence  to  reject  Ho 


14-12-01 


Filename:  undisturbed  ecotypes.xls 
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CEMML  Soil  Data 
Jarvis  Creek  and 
Washington  Range 
Control  Areas 


t-Tests 

Comparing  %  Organic  Matter  Means 


14-12-01 


A  to  D  %  OM 


Scrub  W 

Mixed 

4.98 

5.58 

3.28 

19.10 

3.34 

11.36 

3.38 

24.95 

1.76 

6.63 

3.73 

5.48 

2.73 

12.01 

2.10 

16.52 

0.71 

10.64 

4.78 

7.89 

2.10 

11.54 

3.10 

5.97 

t-Test:  Two-Sample  Assuming  Unequal  Variances 


Scrub  W  Mixed  \ 

Mean 

2.998036 

11.47164 

Variance 

1 .499458 

36.67946 

Observations 

12 

12 

Hypothesized  Mean  Difference 

0 

df 

12 

tStat 

-4.75058 

P(T<=t)  one-tail 

0.000236 

t  Critical  one-tail 

1 .782287 

P(T<=t)  two-tail 

0.000472 

t  Critical  two-tall 

2.178813 

Ho:  Means  are  equal  (difference  =  0) 

Ha:  Means  are  not  equal  (differnce  not  equal  to  0) 
alpha  =  .05 

Conclusion:  Reject  Ho,  means  are  not  equal. 


Filename:  undisturbed  ecotypes.xls 
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CEMML  Soil  Data  Comparing  %  Organic  Matter  Means  for  Craters  and  Plots  14-12-01 

Washington  Impact  Area 
Craters  vs  Plots 
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CEMML  Soil  Data  Comparing  %  Organic  Matter  Means  for  Craters  and  Plots  14-12-01 

Washington  Impact  Area 
and  Control  Area 
Craters  vs  Plots 
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CEMML  Soil  Data  t-Tests 

Crater  Rings  Comparing  %  Organic  Matter  Means  for  Crater  Rings 

Rings  are  Averages  Accross  Crater  7 

Craters 


Comparisons 

Ring  1 

Ring  2 

Ring  3 

Ring  4 

Ring  5 

Ring  1 

X 

X 

X 

X 

Ring  2 

X 

X 

X 

Ring  3 

X 

X 

Ring  4 

X 

Ring  5 

Ring 

Ring  1 

Ring  2 

Ring  3 

Ring  4 

Ring  5 

%  Organic  Matter 

1.00 

0.74 

0.78 

3.46 

4.00 

0.73 

0.86 

1.11 

2.56 

4.88 

0.64 

0.60 

1.02 

2.03 

4.76 

Descriptive  Statistics 

Ring  1 

Ring  2 

Ring  3 

Ring  4 

Ring  5 

Mean 

0.791891 

0.731551 

0.968888 

2.680245 

4.546188 

Standard  Error 

0.107186 

0.073766 

0.0985 

0.417173 

0.276077 

Median 

0.7327 

0.739596 

1.024056 

2.556692 

4.763412 

Mode 

#N/A 

#N/A 

#N/A 

#N/A 

#N/A 

Standard  Deviation 

0.185651 

0.127766 

0.170606 

0.722564 

0.478179 

Sample  Variance 

0.034466 

0.016324 

0.029107 

0.522099 

0.228655 

Kurtosis 

#DiV/0! 

#DlV/0! 

#DIV/0! 

#DIV/0! 

#DIV/0! 

Skewness 

1 .288882 

-0.282238 

-1.302984 

0.74697 

-1 .622374 

Range 

0.356868 

0.255152 

0.32756 

1.429196 

0.87924 

Minimum 

0.643052 

0.599952 

0.777524 

2.027424 

3.997956 

Maximum 

0.99992 

0.855104 

1.105084 

3.45662 

4.877196 

Sum 

2.375672 

2.194652 

2.906664 

8.040736 

13.63856 

Count 

3 

3 

3 

3 

3 

Confidence  Level(95.07o) 

0.461183 

0.317389 

0.42381 

1 .794951 

1.187864 

Summary  of  tests 

Ring  1 

Ring  2 

Ring  3 

Ring  4 

Ring  5  | 

No  evidence  to  reject  Ho: 

Mi— 1 

Means  are  equal  1  1 

■4 — 
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CEMML  Soil  Data  t-Tests 

Crater  Rings  Comparing  %  Organic  Matter  Means  for  Crater  Rings 

Rings  are  Averages  Accross  Crater  7 

Craters 


t-Test:  Two-Sample  Assuming  Unequal  Variances 


Ring  1  Ring  2  \ 

Mean 

0.791891 

0.731551 

Variance 

0.034466 

0.016324 

Observations 

3 

3 

Hypothesized  Mean  Difference 

0 

df 

4 

tStat 

0.46374 

P(T<=t)  one-taii 

0.333473 

t  Critical  one-tail 

2.131846 

P(T<=t)  two-tail 

0.666946 

t  Criticai  two-tail 

2.776451 

Ho:  Means  are  equal  (difference  =  0) 

Ha:  Means  are  not  equal  (differnce  not  equal  to  0) 
alpha  =  .05 

Conclusion:  No  evidence  to  reject  Ho 


t-Test:  Two-Sample  Assuming  Unequal  Variances 


Ring  2  Ring  3  \ 

Mean 

0.731551 

0.968888 

Variance 

0.016324 

0.029107 

Observations 

3 

3 

Hypothesized  Mean  Difference 

0 

df 

4 

tStat 

-1.928644 

P(T<=t)  one-tail 

0.063006 

t  Critical  one-tail 

2.131846 

P(T<=t)  two-tail 

0.126012 

t  Critical  two-tail 

2.776451 

Ho:  Means  are  equal  (difference  =  0) 

Ha:  Means  are  not  equal  (differnce  not  equal  to  0) 
alpha  =  .05 

Conclusion:  No  evidence  to  reject  Ho 


t-Test:  Two-Sample  Assuming  Unequal  Variances 


Ring  3  Ring  4  \ 

Mean 

0.968888 

2.680245 

Variance 

0.029107 

0.522099 

Observations 

3 

3 

Hypothesized  Mean  Difference 

0 

df 

2 

tStat 

-3.992495 

P(T<=t)  one-tail 

0.028694 

t  Criticai  one-tail 

2.919987 

P(T<=t)  two-tail 

0.057388 

t  Critical  two-tail 

4.302656 

Ho:  Means  are  equal  (difference  =  0) 

Ha:  Means  are  not  equal  (differnce  not  equal  to  0) 
alpha  =  .05 

Conclusion:  No  evidence  to  reject  Ho 


t-Test:  Two-Sample  Assuming  Unequal  Variances 


Ring  4  Ring  5  \ 

Mean 

2.680245 

4.546188 

Variance 

0.522099 

0.228655 

Observations 

3 

3 

Hypothesized  Mean  Difference 

0 

df 

3 

tStat 

-3.730009 

P(T<=t)  one-tail 

0.016787 

t  Critical  one-tail 

2.353363 

P(T<=t)  two-tail 

0.033575 

t  Critical  two-tail 

3.182449 

Ho:  Means  are  equal  (difference  =  0) 

Ha:  Means  are  not  equal  (differnce  not  equal  to  0) 
alpha  =  .05 

Conclusion:  Reject  Ho,  means  are  not  equal. 
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CEMMLSoil  Data  t-Tests 

Crater  Rings  Comparing  %  Organic  Matter  Means  for  Crater  Rings 

Rings  are  Averages  Accross  Crater  7 

Craters 


t-Test:  Two-Sample  Assuming  Unequal  Variances 


Ring  1  Ring  3  \ 

Mean 

0.791891 

0.968888 

Variance 

0.034466 

0.029107 

Observations 

3 

3 

Hypothesized  Mean  Difference 

0 

df 

4 

tStat 

-1.215882 

P(T<=t)  one-tail 

0.145439 

t  Critical  one-tall 

2.131846 

P(T<=t)  two-tail 

0.290877 

t  Critical  two-tail 

2.776451 

Ho:  Means  are  equal  (difference  =  0) 

Ha:  Means  are  not  equal  (differnce  not  equal  to  0) 
alpha  =  .05 

Conclusion:  No  evidence  to  reject  Ho 


t-Test:  Two-Sample  Assuming  Unequal  Variances 


Ring  2  Ring  4  | 

Mean 

0.731551 

2.680245 

Variance 

0.016324 

0.522099 

Observations 

3 

3 

Hypothesized  Mean  Difference 

0 

df 

2 

tStat 

-4.599837 

P(T<=t)  one-tail 

0.022078 

t  Critical  one-tail 

2.919987 

P(T<=t)  two-tail 

0.044156 

t  Critical  two-tail 

4.302656 

Ho:  Means  are  equal  (difference  =  0) 

Ha:  Means  are  not  equal  (differnce  not  equal  to  0) 
alpha  =  .05 

Conclusion:  Reject  Ho,  means  are  not  equal. 


t-Test:  Two-Sample  Assuming  Unequal  Variances 


Ring  3  Ring  5  \ 

Mean 

0.968888 

4.546188 

Variance 

0.029107 

0.228655 

Observations 

3 

3 

Hypothesized  Mean  Difference 

0 

df 

3 

tStat 

-12.20412 

P(T<=t)  one-tail 

0.000592 

t  Critical  one-tail 

2.353363 

P(T<=t)  two-tail 

0.001185 

t  Critical  two-tail 

3.182449 

Ho:  Means  are  equal  (difference  =  0) 

Ha:  Means  are  not  equal  (differnce  not  equal  to  0) 
alpha  =  .05 

Conclusion:  Reject  Ho,  means  are  not  equal. 
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CEMML  Soil  Data  t-Tests 

Crater  Rings  Comparing  %  Organic  Matter  Means  for  Crater  Rings 

Rings  are  Averages  Accross  Crater  7 

Craters 


t-Test:  Two-Sample  Assuming  Unequal  Variances 


Ring  1  Ring  4  \ 

Mean 

0.791891 

2.680245 

Variance 

0.034466 

0.522099 

Observations 

3 

3 

Hypothesized  Mean  Difference 

0 

df 

2 

tStat 

-4.384156 

P(T<=t)  one-tail 

0.024145 

t  Critical  one-tail 

2.919987 

P(T<=t)  two-tail 

0.048289 

t  Critical  two-tail 

4.302656 

Ho:  Means  are  equal  (difference  =  0) 

Ha:  Means  are  not  equal  (differnce  not  equal  to  0) 
alpha  =  .05 

Conclusion:  Reject  Ho,  means  are  not  equal. 


t-Test:  Two-Sample  Assuming  Unequal  Variances 


Ring  2  Ring  5  j 

Mean 

0.731551 

4.546188 

Variance 

0.016324 

0.228655 

Observations 

3 

3 

Hypothesized  Mean  Difference 

0 

df 

2 

tStat 

-13.34901 

P(T<=t)  one-tail 

0.002782 

t  Critical  one-tail 

2.919987 

P(T<=t)  two-tail 

0.005565 

t  Critical  two-tail 

4.302656 

Ho:  Means  are  equal  (difference  =  0) 

Ha:  Means  are  not  equal  (differnce  not  equal  to  0) 
alpha  =  .05 

Conclusion:  Reject  Ho,  means  are  not  equal. 


t-Test:  Two-Sample  Assuming  Unequal  Variances 


Ring  1  Ring  5\ 

Mean 

0.791891 

4.546188 

Variance 

0.034466 

0.228655 

Observations 

3 

3 

Hypothesized  Mean  Difference 

0 

df 

3 

tStat 

-12.67684 

P(T<=t)  one-tail 

0.000529 

t  Critical  one-tail 

2.353363 

P(T<=t)  two-tall 

0.001059 

t  Critical  two-tail 

3.182449 

Ho:  Means  are  equal  (difference  =  0) 

Ha:  Means  are  not  equal  (differnce  not  equal  to  0) 
alpha  =  .05 

Conclusion:  Reject  Ho,  means  are  not  equal. 
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%  Organic  Matter 


Descriptive  Statistics 


Mean 


Standard  Error 


Median 


Mode 


Standard  Deviation 


Sample  Variance 


Kurtosis 


Skewness 


Range 


Minimum 


Maximum 


Count 


Confidence  Level(95.0%) 


Ring  2 


X 


ffiBEI 

liEEi 

SJEQI 


Ring  3 


X 


X 


Ring  4 


X 


X 


X 


Ring  5 


X 


X 


X 


X 


Ring  1 

Ring  2 

Ring  3 

Ring  4 

Ring  5 

1.57 

0.91 

0.60 

0.93 

1.18 

1.43 

1.52 

1.55 

1.23 

1.51 

1.43 

1.31 

0.56 

1.04 

1.18 

Ring  1 

Ring  2 

Ring  3 

Ring  4 

Ring  5 

1.475169 

1.247601 

0.903951 

1 .063708 

1.288403 

0.045136 

0.181051 

0.323183 

0.087432 

0.108353 

1,432644 

1.311964 

0.601676 

1.039572 

1.184388 

#N/A 

#N/A 

#N/A 

#N/A 

#N/A 

0.078178 

0.31359 

0.55977 

0.151437 

0.187673 

0.006112 

0.098338 

0.313343 

0.022933 

0.035221 

#DlV/0! 

#DlV/0! 

#DIV/0! 

#DIV/0! 

#DIV/0! 

1 .723526 

-0.884695 

1.721411 

0.698989 

1.727941 

0.13792 

0.617192 

0.989576 

0.299976 

0.329284 

1.427472 

0.906824 

0.5603 

0.925788 

1.175768 

1.565392 

1.524016 

1 .549876 

1 .225764 

1.505052 

4.425508 

3.742804 

2.711852 

3.191124 

3.865208 

3 

3 

3 

3 

3 

0.194205 

0.779 

1.390547 

0.376192 

0.466207 

Ring  1  | 

1  Ring  2  | 

I  Ring  3 

[  Ring  4 

1  Ring  5  | 

14-12-01 


CEMML  Soil  Data  t-Tests 

Crater  Rings  Comparing  %  Organic  Matter  Means  for  Crater  Rings 

Rings  are  Averages  Accross  Crater  1 7 

Craters 


t-Test:  Two-Sample  Assuming  Unequal  Variances 


Ring  1  Ring  2  \ 

Mean 

1.475169 

1.247601 

Variance 

0.006112 

0.098338 

Observations 

3 

3 

Hypothesized  Mean  Difference 

0 

df 

2 

tStat 

1.219599 

P(T<=t)  one-taii 

0.173461 

t  Critical  one-tail 

2.919987 

P(T<=t)  two-tail 

0.346922 

t  Criticai  two-taii 

4.302656 

Ho:  Means  are  equal  (difference  =  0) 

Ha:  Means  are  not  equal  (differnce  not  equal  to  0) 
alpha  =  .05 

Conclusion:  No  evidence  to  reject  Ho 


t-Test:  Two-Sample  Assuming  Unequal  Variances 


1  Ring  2  Ring  3  | 

Mean 

1.247601 

0.903951 

Variance 

0.098338 

0.313343 

Observations 

3 

3 

Hypothesized  Mean  Difference 

0 

df 

3 

tStat 

0.927678 

P(T<=t)  one-taii 

0.211 

t  Critical  one-tail 

2.353363 

P(T<=t)  two-tail 

0.422001 

t  Critical  two-tail 

3.182449 

Ho:  Means  are  equal  (difference  =  0) 

Ha:  Means  are  not  equal  (differnce  not  equal  to  0) 
alpha  =  .05 

Conclusion:  No  evidence  to  reject  Ho 


t-Test:  Two-Sample  Assuming  Unequal  Variances 


Ring  3  Ring  4  \ 

Mean 

0.903951 

1 .063708 

Variance 

0.313343 

0.022933 

Observations 

3 

3 

Hypothesized  Mean  Difference 

0 

df 

2 

tStat 

-0.47717 

P(T<=t)  one-tail 

0.340149 

t  Critical  one-tail 

2.919987 

P(T<=t)  two-tail 

0.680298 

t  Critical  two-tail 

4.302656 

Ho:  Means  are  equal  (difference  =  0) 

Ha:  Means  are  not  equal  (differnce  not  equal  to  0) 
alpha  =  .05 

Conclusion:  No  evidence  to  reject  Ho 


t-Test:  Two-Sample  Assuming  Unequal  Variances 


Ring  4  Ring  5^ 

Mean 

1.063708 

1 .288403 

Variance 

0.022933 

0.035221 

Observations 

3 

3 

Hypothesized  Mean  Difference 

0 

df 

4 

tStat 

-1.613843 

P(T<=t)  one-tail 

0.09093 

t  Criticai  one-taii 

2.131846 

P(T<=t)  two-tail 

0.18186 

t  Critical  two-tail 

2.776451 

Ho:  Means  are  equal  (difference  =  0) 

Ha:  Means  are  not  equal  (differnce  not  equal  to  0) 
alpha  =  .05 

Conclusion:  No  evidence  to  reject  Ho 
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CEMML  Soil  Data  t-Tests 

Crater  Rings  Comparing  %  Organic  Matter  Means  for  Crater  Rings 

Rings  are  Averages  Accross  Crater  1 7 

Craters 


t-Test:  Two-Sample  Assuming  Unequal  Variances 


Ring  1  Ring  3  \ 

Mean 

1.475169 

0.903951 

Variance 

0.006112 

0.313343 

Observations 

3 

3 

Hypothesized  Mean  Difference 

0 

df 

2 

tStat 

1 .750486 

P(T<=t)  one-tail 

0.111068 

t  Critical  one-tail 

2.919987 

P(T<=t)  two-tail 

0.222137 

t  Critical  two-tail 

4.302656 

Ho:  Means  are  equal  (difference  =  0) 

Ha:  Means  are  not  equal  (differnce  not  equal  to  0) 
alpha  =  .05 

Conclusion:  No  evidence  to  reject  Ho 


t-Test:  Two-Sample  Assuming  Unequal  Variances 


Ring  2  Ring  4  \ 

Mean 

1 .247601 

1.063708 

Variance 

0.098338 

0.022933 

Observations 

3 

3 

Hypothesized  Mean  Difference 

0 

df 

3 

tStat 

0.914633 

P(T<=t)  one-tail 

0.213914 

t  Critical  one-tall 

2.353363 

P(T<=t)  two-tail 

0.427828 

t  Critical  two-tall 

3.182449 

Ho:  Means  are  equal  (difference  =  0) 

Ha:  Means  are  not  equal  (differnce  not  equal  to  0) 
alpha  =  .05 

Conclusion:  No  evidence  to  reject  Ho 


t-Test:  Two-Sample  Assuming  Unequal  Variances 


Ring  3  Ring  5  \ 

Mean 

0,903951 

1.288403 

Variance 

0.313343 

0.035221 

Observations 

3 

3 

Hypothesized  Mean  Difference 

0 

df 

2 

tStat 

-1.127876 

P(T<=t)  one-tail 

0.188242 

t  Critical  one-tail 

2.919987 

P(T<=t)  two-tail 

0.376484 

t  Critical  two-tail 

4.302656 

Ho:  Means  are  equal  (difference  =  0) 

Ha:  Means  are  not  equal  (differnce  not  equal  to  0) 
alpha  =  .05 

Conclusion:  No  evidence  to  reject  Ho 
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CEMML  Soil  Data  t-Tests 

Crater  Rings  Comparing  %  Organic  Matter  Means  for  Crater  Rings 

Rings  are  Averages  Accross  Crater  1 7 

Craters 


t-Test:  Two-Sample  Assuming  Unequal  Variances 


Ring  1  Ring  4  \ 

Mean 

1.475169 

1,063708 

Variance 

0.006112 

0.022933 

Observations 

3 

3 

Hypothesized  Mean  Difference 

0 

df 

3 

tStat 

4.181704 

P(T<=t)  one-tail 

0.012459 

t  Critical  one-tail 

2.353363 

P(T<=t)  two-tail 

0.024918 

t  Critical  two-tail 

3.182449 

Ho:  Means  are  equal  (difference  =  0) 

Ha:  Means  are  not  equal  (differnce  not  equal  to  0) 
alpha  =  .05 

Conclusion:  Reject  Ho,  means  are  not  equal. 


t-Test:  Two-Sample  Assuming  Unequal  Variances 


Ring  2  Ring  5  \ 

Mean 

1.247601 

1 .288403 

Variance 

0.098338 

0.035221 

Observations 

3 

3 

Hypothesized  Mean  Difference 

0 

df 

3 

tStat 

-0.193374 

P(T<=t)  one-tail 

0.429509 

t  Critical  one-tall 

2.353363 

P(T<=t)  two-tail 

0.859018 

t  Critical  two-tail 

3.182449 

Ho:  Means  are  equal  (difference  =  0) 

Ha:  Means  are  not  equal  (differnce  not  equal  to  0) 
alpha  =  .05 

Conclusion:  No  evidence  to  reject  Ho 


t-Test:  Two-Sample  Assuming  Unequal  Variances 


Ring  1  Ring  5  \ 

Mean 

1.475169 

1.288403 

Variance 

0.006112 

0.035221 

Observations 

3 

3 

Hypothesized  Mean  Difference 

0 

df 

3 

tStat 

1.59115 

P(T<=t)  one-tail 

0.104904 

t  Critical  one-tall 

2.353363 

P(T<=t)  two-tail 

0.209808 

t  Critical  two-tail 

3.182449 

Ho:  Means  are  equal  (difference  =  0) 

Ha:  Means  are  not  equal  (differnce  not  equal  to  0) 
alpha  =  .05 

Conclusion:  No  evidence  to  reject  Ho 
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